AD-R158  993 
UNCLASSIFIED 


COMPUTATION  OF  STRESS  INTENSITY  FACTORS(U)  HASHINQTON 
UNIV  ST  LOUIS  NO  CENTER  FOR  COMPUTATIONAL  MECHANICS 
B  A  SZABO  ET  AL.  OCT  B3  HU/CCH-82/3  N00014-81-K-0S25 

F/G  20/11 


1/2 


NL 


**■ 


WASHINGTON 
UNIVERSITY 
IN  ST  LOUIS 


REPORT  WU/CCM-83/3 


00 

O) 

o> 

CO 

in 

< 

i 

O 

< 


COMPUTATION  OF 
STRESS  INTENSITY 
FACTORS 


Barna  A.  Szabo 
and 

Dionysios  Vasilopoulos 


OCTOBER,  1983 


Prepared  for: 

Office  of  Naval  Research 
Arlington,  VA  2221 7 


S 


DTIC 

ELECTE 
SEP  1  1  1985 


r 


CENTER  FOR 
COMPUTATIONAL 
MECHANICS 

WASHINGTON  UNIVERSITY 
CAMPUS  BOX  1129 
ST.  LOUIS.  MO  63130 


|  TaU  docu*o«U  boa  b—e  uppiot»d  | 
fat  public  r«i«oa*  end  aerial  Hi 
distribution  fa  uliaML 


ABSTRACT 


I 

t. 

^  The  performance  characteristics  of  the  generalized 
influence  function  method  for  the  approximate  computation  of 
the  amplitudes  of  the  eigenfunctions  of  the  equations  of 
plane  elasticity  in  the  vicinity  of  sharp  reentrant  corners 
were  evaluated.  The  eigenfunctions  satisfy  the  equations  of 
equilibrium,  compatibility  and  stress-strain  laws  and  the 
free-free  boundary  conditions  at  reentrant  corners.  The 
amplitudes  of  the  eigenfunctions  are  called  the  generalized 
stress  intensity  factors. 


It  is  concluded  that  the  generalized  stress  intensity 
factors  can  be  computed  to  within  one  percent  relative  error 
with  small  computational  effort.  Therefore  the  essential 
characteristics  of  the  elastic  stress  field  in  the  neighbor¬ 
hood  of  reentrant  corners  can  be  determined  with  great 
precision.  This  computational  technology  is  essential  for 


the  development  of  theories  of  crack  initiation  in  metals 
and  composites.  -  —  ,  /  .  ■  /  . 
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TREATMENT  OF  GEOMETRIC  SINGULARITIES  WITH 
THE  p- VERS ION  OF  THE  FINITE  ELEMENT  METHOD 

1 .  INTRODUCTION 

In  the  displacement  formulation  of  the  finite  element 
method  the  solution  minimizes  the  strain  energy  of  the 
error  for  the  given  finite  element  mesh  and  polynomial 
degree  of  elements.  It  has  been  shown  that  this  is  closely 
related  to  minimizing  the  root-mean-square  error  in  stress 
[1].  In  engineering  computations,  however,  the  strain  energy 
and  the  root-mean- square  error  in  stress  are  not  the  quanti¬ 
ties  of  primary  interest.  The  goal  of  computation  is  usually 
to  estimate  within  a  reasonably  small  margin  of  error,  some 
functional  of  the  displacement  field,  such  as  stress  at  a 
point,  stress  intensity  factors  etc.  These  quantities  are 
usually  computed  from  the  finite  element  solution  directly, 
without  taking  into  account  that  the  essential  character  of 
the  solution  is  known  a  priori.  For  example,  stresses  at  a 
point  within  an  element  are  computed  as  some  linear 

*The  numbers  in  brackets  in  the  text  indicate  references 
in  the  Bibliography. 


-2- 


combination  of  the  derivatives  of  the  shape  functions  at 
the  point.  Stress  intensity  factors  for  crack  problems 
are  computed  from  the  strain  energy  release  rate,  which  is 
obtained  by  computing  the  strain  energy  for  two  bodies  with 
slightly  different  crack  sizes.  Although  the  method  per¬ 
forms  well  when  used  in  conjunction  with  the  p- version  of 
the  finite  element  method  [2],  the  stress  intensity  factors 
for  mixed  mode  problems  cannot  be  obtained  separately  by 
this  method. 

In  three  related  papers  Babuska  and  Miller  introduced 
new  techniques  that  permit  the  extraction  of  various 
functionals  of  the  displacements,  such  as  pointwise 
displacements,  stresses  and  stress  intensity  factors,  with 
greater  reliability  and  accuracy  than  was  previously 
possible  [3,  4,  5].  The  main  idea  of  the  new  extraction 
techniques  is  that  not  only  the  finite  element  solution  J 

but  also  the  essential  characteristics  of  the  exact 
solution  are  utilized  in  making  the  computations.  This  *j 

requires  a  modest  amount  of  extra  effort  but  a  great  deal  ^ 

is  gained  in  accuracy  and  reliability.  One  of  the  extrac-  1 

tion  techniques  is  the  generalized  influence  function 
method  for  the  computation  of  stress  intensity  factors. 

By  stress  intensity  factors  we  mean  not  only  the  stress 
intensity  factors  defined  in  linear  elastic  fracture 


mechanics  but  also  the  analogous  quantities  associated 
with  any  sharp  reentrant  corner. 

Stress  singularities  in  the  solution  of  elasticity 
problems  can  be  caused  by  loading,  sudden  changes  in  the 
boundary  conditions  or  material  properties  and  reentrant 
corners.  A  typical  corner  detail  is  shown  in  figure  1.1. 
This  investigation  is  concerned  with  stress  singularities 
at  reentrant  corners  with  free-free  boundary  conditions 
on  the  sides  of  the  angle. 

In  the  neighborhood  of  the  corner  the  solution  vector 
can  be  written  in  terms  of  polar  coordinates  centered  on 
the  corner  in  the  following  form: 

<  . 

u  -  Z  Ki  r  1  F i  ( 9 )  +  G  (r ,  8)  (1.1) 

where  u  ■  (u^,  Uj) ,  F  and  G  are  functions  that  are  smoother 
than  r  1  F^(8)  and  the  <^'s  are  positive  numbers.  The 
amplitudes  are  the  generalized  stress  intensity  factors. 
The  eigenvalues  depend  on  the  angle  6,  the  boundary 
conditions  imposed  on  the  two  sides  of  the  angle,  and 
Poisson's  ratio.  The  eigenfunction  F^(9)  depends  on  the 
angle  2,  the  corresponding  eigenvalue  and  the  elastic 
constants.  A  typical  stress  component  is  of  the  form: 


-  Z  K.  r 

i 

i 


f(r)  +  g (r ,  9) 


(1.2) 
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It  can  be  seen  that  for  values  of  < .  <  1  tha  stresses 

1 

become  infinite  at  the  corner. 

A  special  case  is  the  case  of  a  crack  (angle  8 
equal  to  360  degrees) .  The  solution  of  crack  problems 
and  the  study  of  the  conditions  under  which  cracks 
propogate  is  the  subject  of  linear  elastic  fracture 
mechanics.  Its  applicability  rests  on  the  satisfaction 
of  a  criterion  known  as  small  scale  yielding.  This  means 
that  any  nonlinearities  are  confined  to  a  region  which  is 
small  in  comparison  with  the  size  of  the  body  and  is  com¬ 
pletely  surrounded  by  a  region  in  which  the  solution  of  the 
elasticity  problem  is  an  adequate  representation  of  the  real 
response  of  the  material.  The  elastic  stress  field  and 
strain  field  in  the  neighborhood  of  the  singularity  can  be 

written  respectively  as: 

K  ** 

c..  ■  -  a. . (e)  +  higher  order  terms  (1.3) 


V 

£..  »  ■  ■  e..(8)  +  higher  order  terms  (1.4) 

^  /Jnr 

where  K  is  the  stress  intensity  factor  and  0^(8), 

£^ (8)  are  smooth  functions  of  8.  When  the  restrictions 
stated  in  the  preceding  paragraph  are  met,  these  fields  are 
the  same  for  any  crack,  regardless  of  the  overall  geometry 
and  applied  loading.  The  stress  and  deformation  fields 
ahead  of  the  crack  tip  for  a  purely  symmetric  or 


antisymmetric  crack  configuration  are  then  characterized 
by  a  single  parameter,  the  stress  intensity  factor  K. 

Linear  elastic  fracture  mechanics  is  based  on  the  obser¬ 
vation  that  crack  growth  is  controlled  by  the  parameter 
K.  A  pre-cracked  specimen  in  which  the  elastic  field  is 
known  is  tested  under  monotonic  loading  conditions  until 
the  first  occurence  of  crack  propagation  is  observed. 

At  this  point  the  stress  intensity  factor  is  said  to  have 
a  critical  value  which  is  usually  denoted  by  Kc.  This 
phenomenological  approach  ignores  the  micromechanisms  of 
void  nucleation  ahead  of  the  crack  tip  and  the  way  in 
which  these  voids  are  joined  to  increase  the  crack  length, 
yet  it  has  been  highly  successful  in  giving  reliable 
answers  to  important  practical  questions  such  as  the 
expected  life  of  structures  under  existing  flaws,  the 
maximum  flaw  size  that  quality  control  is  allowed  to  miss 
and  maximum  allowable  time  intervals  between  inspections . 
The  critical  stress  intensity  factor  KJC  in  mode  I  under 
plane  strain  conditions  is  termed  fracture  toughness  and 
is  accepted  as  the  material  constant  that  characterizes 
material  resistance  to  fracture. 

This  approach  has  been  extended  to  the  nonlinear 
regime.  The  J-integral  is  the  intensity  of  the  singular 
nonlinear  elastic  strain  field  in  the  neighborhood  of  the 
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of  the  crack  tip  (HRR  field) .  In  this  sense  J  c  is 
analogous  to  KJC  but  of  course  the  small  scale  yielding 
condition  is  now  replaced  by  the  " J-dominance"  condition, 
which  means  that  the  area  ahead  of  the  tip,  in  which 
intense  nonlinearities  deviating  from  deformation  theory 
of  plasticity  take  place,  is  relatively  small  and  is 
completely  surrounded  by  a  region  in  which  the  singular 
nonlinear  elastic  field  represents  adequately  the  real 
material  response  [6,  7,  8,  9,  10]. 

Despite  its  great  success,  linear  elastic  fracture 
mechanics  is  difficult  to  apply  to  structures  of  compli¬ 
cated  geometry.  The  presence  of  a  large  number  of  points 
where  cracks  are  apt  to  occur  would  require  a  large 
number  of  analyses  to  be  performed  with  a  hypothetical 
crack  at  each  of  these  points .  On  the  other  hand ,  the 
designer's  intent  is  that  the  structure  should  spend  most 
of  its  life  in  the  crack  initiation  stage  rather  than  in 
the  crack  propagation  stage.  The  ability  to  formulate 
crack  initiation  criteria  based  on  linear  elasticity 
could  open  the  way  to  more  rational  design  procedures. 

Although  the  singular  nature  of  elasticity  solutions 
in  the  neighborhood  of  reentrant  corners  was  reported 
in  the  engineering  literature  as  early  as  1933  [ll], 


•  A 
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and  the  1952  paper  by  M.  L.  Williams  [12]  gave  the  strength 
of  these  singularities  for  the  various  angles,  no  attempt 
was  previously  made  to  assign  any  physical  significance  to 
the  amplitudes  of  the  stress  singular  terms  (except  in 
linear  elastic  fracture  mechanics)  and  no  numerical  method 
existed  for  the  computation  of  these  quantities.  Based  on 
the  work  presented  herein,  these  amplitudes  can  now  be 
computed  with  levels  of  precision  normally  expected  in 
engineering  computations  at  the  expenditure  of  a  relatively 
modest  computational  effort. 

Our  ability  to  compute  the  amplitude  of  all  terms  of 
the  asympotic  expansion  in  the  neighborhood  of  a  reentrant 
corner  of  any  size  offers  new  possibilities  in  the  area  of 
failure  initiation.  When  sufficient  number  of  terms  are 
used  in  the  expansion  and  the  amplitudes  are  accurately 
computed,  then  all  stress  field  parameters  are  known  in 
the  neighborhood  of  corner  points,  therefore  various 
hypotheses  concerning  relationships  between  elastic  stress 
field  parameters  and  failure  initiation  can  be  tested.  in 
the  absence  of  proper  extracton  methods,  uncertainties  in 
numerically  computed  stress  field  parameters  render  the 
formulation,  testing  and  applications  of  such  hypotheses 
very  tenuous . 

A  possible  hypothesis  for  example  is  that  the 
generalized  stress  intensity  factors  are  responsible  for 
crack  initiation  in  reentrant  corners  in  the  same  way  that 
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Let  us  consider  the  case  of  mode  I  and  rewrite 
equations  (2.16)  and  (2.17)  in  the  following  form: 


u  =  K  r1/2  GO) 


v  =.  K  r1//2  H  (9) 


(2.20) 

(2.21) 


Their  derivatives  can  be  expressed  in  the  form: 


4=  -  Kj  r"1/2  G1(0);  ||  -  Kj  r-1/2  G2  C6) 


3v 


(2.22) 


3  v 
3x 


Kj  r“1/2  H1(6) 


3  v 
3y 


Kj  r“1/2  H2(0) 


(2.23) 


We  choose  the  auxiliary  functions  to  be  of  the  following 


form:  ’  .. 

•1 

«  -  r_1/2  <M0)  (2.25) 

*  r"1/2  Y  (0)  .  (2.26) 

Their  derivatives  can  be  expressed  as:  - 

§|  =  r_3/2  $1(6) ,  H  «  r"3/2  *2(0)  (2.27) 

y  j 

||  *  r_3/2  4'1(0)  ,  ||  =  r'3/2  H^O)  .  (2.28)  - 


4 
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the  displacements  u  and  v  is  known  [13,  14]  and  is  given 
in  the  case  of  plane  strain  by: 


Gu  =  Kj  (2tt)“1/2  r1/2  cos  |  (2-2v-cos2  |  ) 

Gv  =  KI  (2tt)"1/2  r1/2  sin  |  (2-2v-cos2  |  ) 
for  mode  I  (opening  mode),  and  by: 

Gu  «  KXI  (2tt)“1/2  r1/2  sin  |  (2-2v+cos2  |  ) 

Gv  =  Kjj  (2,)-1/2  r1//>2  cos  (2v-cos2  ) 


(2.16) 

(2.17) 


(2.18) 

(2.19) 


for  mode  XI  (sliding  mode) . 

These  are  only  the  first  terms  in  the  asymptotic 
expansions.  The  second  terms  are  of  the  order  r.  The 
third  terms  are  of  the  order  r2//2  and  so  on.  The 
asymptotic  expansions  also  contain  the  corresponding 
negative  powers  of  r  (-1/2,  -1,  -3/2  etc.),  but  these 
terms  make  no  sense  physically  since  they  would  imply 
infinite  displacements  at  the  crack  tip,  so  they  are 
not  considered.  Nevertheless,  these  eigenfunctions 
corresponding  to  the  negative  eigenvalues  will  be  very 
useful  for  our  formulation:  the  auxiliary  functions 
0  and  ^  are  chosen  to  be  precisely  these  eigenfunctions. 
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In  view  of  (2.3)  and  (2.4)  we  can  write  (2.12)  and 
(2.13)  as: 

I5  *  X  <j)  (2.14) 

Ig  *  Y  \p.  (2.15) 

On  that  part  of  the  boundary  where  surface  tractions  are 
specified,  the  quantities  Ig  and  Ig  are  known  and  their 
integrals  are  computed  numerically .  On  that  part  of 
the  boundary  where  displacements  are  specified,  the  $ 
and  \p  functions  vanish.  The  reason  that  the  functions 
$  and  tp  are  chosen  to  vanish  there  is  that  we  do  not 
wish  the  extraction  expression  to  contain  a  contour 
integral  of  the  derivatives  of  displacement  components. 

The  extraction  method  has  a  rate  of  convergence  equal 
to  the  rate  of  convergence  in  energy  if  the  integrals 
used  in  the  extraction  represent  only  energy  expressions . 
This  means  that  area  integrals  may  contain  up  to  first 
order  derivatives  of  the  displacements ,  whereas  contour 
integrals  may  contain  only  the  displacements  themselves. 

It  is  observed  that  expressions  Ig  and  I4  fall  into  this 
category  and  so  they  do  not  reduce  the  rate  of  convergence. 

In  order  to  evaluate  the  contour  integrals  on  the 
circular  arc  r  we  proceed  as  follows:  in  the  neighborhood 
of  the  crack  where  the  circular  arc  is  located,  the  form  of 


20 


+  G(  If  nx  +  If  ny>]v 


(2.11) 


W*‘li*S?>»*  +  e‘ii».  +  i7V  + 


+  G(  ft  nx  +  §7  V]0 


(2.12) 


h/  3v  i  *  rf  3v  «  ,  Sv  .  . 

[X<  5S  +  57  1  +  G(  5J  nx  +  57  ny>  + 


+  G'-57ny  +  57V^ 


(2.13) 


In  (2.8)  we  recognize  the  expression  in  the  bracket 
as  the  force  in  the  x  direction  in  terms  of  displace¬ 
ments  (<J>,  iJj)  .  Similarly  in  (2.9)  the  bracket  represents 
the  force  in  the  y  direction  corresponding  to  those 
displacements  (<p,  ip)  .  If  the  auxiliary  functions  4>  and 
are  chosen  to  satisfy  the  equations  of  equilibrium, 
then  I^  and  I 2  vanish. 
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•5 


£  (X+G)  |f  vny  ♦  X  £  unx  +  G  g  nxv  + 


9  tii 


3<|i 


r+r 


+  G  If  nyv  +  G  H  uny]ds 


^  [(X+G)  |s  *nv  +  X  |H  *„v  +  G  ||  V  + 


r+r 


+  G  |f  ny*  +  G  |f  nx»]ds 


=  0 


(2.6) 


Adding  (2.5)  and  (2.6)  we  obtain: 


U 


<I2+I2)dA  -  ( 

r+r. 


(I3+I4)ds 


•5 


(Ig+Ig) ds=0 


r+r. 


(2.7) 


where : 


+  +Gt^  +  0,]U 


(2.8) 


*2-  E'**G)  £  <  !£  +  M>  +G(  0  +  -&,]v  <2-9) 


I, -[X(M  +  If 1  n*  +  s(!£n-  +  l£n 


3x 


3y  y 


•1 
*  « 

*  mf 


N 


3 


+  G  ( 


n 


+ 


3ti 

3x 


n  ]u 

y 


(2.10) 


and,  finally: 


dX  dX  9y  dy 


+  G  Sill 


H  +  If  >*"*  +  G  Sf  *"x  +  G  If  *ny  + 


9y  T  y 


+  G  fx  ^“x  +  G  Tx  4>ny^dS 


S“d“*  4 


V  +  Gil*u  +  Gi!*u  + 


+  G  fx^  u  +  G  ^x3y.v-'dA 


[x  H  u  nx  +  x  ft  v  ny  +  5  H  u  nx  + 


+Gf^uny+Gfxunx+Gfyv  "x^5 


tx  H  ♦  nx  +  X  w  *  nx  +  G  Si  *  nx  + 


6  ♦  n. .  +  G  u  M.  +  G  I?  ♦  n, ]ds 
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3  /  3u  3v  .  r.  32v  32v  >  . 

(A+G)  17  (  *x  +  I?  }  +  G{  +  >  -  °' 


(2.2) 


The  stress  boundary  conditions  are: 


rr  .  ,  3u  ,  3v  .  J  ,  3u  ^  3u  . 
X  =  A(  3x  3y  )nx  G(  3x  nx  3y  ny) 


+  s(  i  n*  +  &  y 


(2.3) 


.  .  3u  3v  .  3v  3v  . 

A(  lx  +  37  )ny  +  G<  3x  nx  +  37  ny}  + 


+  G(  If  nx  +  If  ny) 


(2.4) 


where  X,  g  are  Lame's  constants  and  n  *  (nx,  ny) 
is  the  outward  normal  to  the  surface.  X  and  Y  are 
the  vector  components  of  specified  tractions  in  the  x  and 
y  directions.  Let  us  now  choose  two  functions  $  and  ip , 
the  properties  of  which  will  be  discussed  later  and  let 
us  multiply  equation  (2.1)  by  $  and  integrate  over  the 
domain.  Applying  the  Gauss  theorem  twice  we  obtain: 


Ja  $1*  £<$♦$? 


)■  +  G( ,  + 


3x  3y 


2  2 
.  r  3  u  .  -  3  v 

rr  3x3y 


J  <}>dA 
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2.  THE  GENERALIZED  INFLUENCE  FUNCTION  METHOD 
IN  LINEAR  ELASTIC  FRACTURE  MECHANICS 

We  consider  a  two  dimensional  body  containing  a  crack 
as  shown  in  figure  (2.1).  The  boundary  consists  of  piece- 
wise  smooth  curves  and  there  is  no  re-entrant  corner  other 
than  the  one  at  the  crack  tip.  Extension  to  cases  where 
the  domain  contains  more  than  one  geometric  singularity 
will  not  introduce  any  additional  difficulties.  This  will 
be  shown  later. 

We  now  remove  from  the  domain  a  small  disk  of  radius 
r»e  centered  on  the  crack  tip  as  shown  in  figure  (2.2) . 

The  boundary  now  consists  of  two  parts:  the  circular  arc 
r£  and  the  rest  of  the  original  boundary  r.  The  removal 
of  the  singularity  enables  us  to  perform  integration  by 
parts.  The  limit  is  then  taken  as  r  tends  to  zero. 

Let  us  choose  a  Cartesian  coordinate  system  ( x,y ) 
centered  on  the  crack  tip  with  the  x-axis  in  the  direction 
of  the  crack ,  and  the  corresponding  polar  system  (r,6)  as 
shown  in  figure  (2.3),  and  write  the  equations  of  equili¬ 
brium  in  terms  of  displacements  (u,v)  in  the  Cartesian 
system: 

U+G)  (  +  U  }  +G(^f+^r)  =  0  (2.1) 

3X  3X  3y  3x2  3y2 
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In  chapter  7  it  is  concluded  that  the  method  is 


separate  the  stress  intensity  factors  for  inodes  I  and  II 
is  demonstrated. 

In  chapter  3  explicit  expressions  ai re  derived  for 
the  asymptotic  expansions  of  the  displacements  in  the 
neighborhood  of  reentrant  corners  of  any  size.  The  real 
eigenvalues  corresponding  to  modes  I  and  II  for  various 
angles  are  tabulated.  The  displacement  eigenfunctions 
are  also  derived  for  the  case  of  complex  eigenvalues. 

In  chapter  4  the  generalized  influence  function 
method  is  implemented  for  the  case  of  a  reentrant  corner 
of  arbitrary  size.  The  method  for  obtaining  the  amplitude 
of  any  term  in  the  expansion  is  demonstrated.  Special 
consideration  is  given  to  the  case  of  complex  eigenvalues .  ■ 

In  chapter  5  the  implementation  is  tested  against  a 
model  problem  for  which  the  exact  solution  is  known.  This 
allows  rigorous  convergence  study  to  be  performed.  The 
theoretically  predicted  rates  of  convergence  for  various 
reentrant  corners  are  verified  numerically. 

In  chapter  6  tests  performed  on  double  edge  notched 
epoxy  specimens  are  discussed.  This  test  data  shows  that 
a  monotonic  relationship  exists  between  the  stress  intensity 
factor  and  failure  initiation  for  a  wide  range  of  solid 
angles. 


function  method  for  the  extraction  of  stress  intensity 
factors  in  plane  elasticity - 

In  the  case  of  linear  elastic  fracture  mechanics  the 
generalized  influence  function  method  yields  separately 
the  stress  intensity  factors  corresponding  to  the  symmetric 
mode  of  deformation  (mode  I)  and  to  the  antisymmetric  mode 
(mode  II) .  Our  ability  to  separate  of  the  two  modes  may 
prove  to  be  important  in  linear  elastic  fracture  mechanics, 
where  only  the  combination  of  the  two  modes  could  be  com¬ 
puted  previously  with  reasonable  accuracy  and  the  contribu¬ 
tion  of  mode  II  has  not  been  well  understood. 

The  method  has  a  rate  of  convergence  equal  to  the 
rate  of  convergence  in  energy,  which  is  twice  the  rate  of 
convergence  in  energy  norm  [4].  Our  ability  to  achieve  con 
vergence  in  practical  computations  can  be  used  as  a  tool 
for  error  estimation.  Knowing  the  theoretical  rate  of 
convergence,  we  can  obtain  an  estimate  of  the  exact  value 
of  the  stress  intensity  factor  (also  strain  energy,  root- 
mean  square  stress  measure  and  other  functionals)  accurate 
to  within  one  percent  relative  error. 

In  chapter  2  the  method  is  implemented  for  crack 
problems.  The  rationale  of  the  method  is  explained  and 
the  extraction  functions  are  derived.  The  ability  to 


the  stress  intensity  factor  in  linear  elastic  fracture 
mechanics  is  responsible  for  crack  propagation.  Again, 
some  small  scale  yielding  criterion  must  be  satisfied. 

If  experiments  verify  this  hypothesis  and  establish  its 
limits  of  applicability,  then  we  have  at  our  disposal  a 
very  simple  and  powerful  method  for  the  design  of  structures 
with  geometric  singularities .  The  methodology  of  linear 
elastic  fracture  mechanics  can  be  generalized  in  a 
straightforward  manner.  For  a  given  angle  and  mode  of 
loading  the  singular  elastic  field  in  the  neighborhood  of 
a  reentrant  corner  is  of  the  same  form  regardless  of  the 
overall  configuration  and  load  distribution.  Any  change  in 
the  loading  or  the  boundary  conditions  affects  material 
behavior  at  the  tip  of  the  notch  only  through  the 
generalized  stress  intensity  factors  which  are  the  ampli¬ 
tudes  of  the  terms  in  the  asymptotic  expansion  of  the 
linear  elastic  solution  in  the  neighborhood  of  reentrant 
corners  of  arbitrary  size.  For  large  notch  angles  the 
higher  order  terms  may  also  be  important  as  crack 
initiation  parameters.  We  can  then  argue  that  under 
conditions  of  small  scale  yielding  the  generalized  stress 
intensity  factors  can  be  used  for  predicting  failure 
initiation  events. 

The  scope  of  this  report  is  the  implementation, 
application  and  evaluation  of  the  generalized  influence 


The  integrals  of  1^,  I 4 /  I 5  and  Ig  along  the  circular 
arc  can  now  be  computed.  A  typical  term  for  the 
integral  of  I 3  will  be: 

j  H  nx  u  ds 

r 

e 

substituting 

II  -  r“3/2  *1(e),  nn  -  -cose,  u  *  <I  r1/2  6(0), 

ds  =  rd6 


we  obtain 


I 


5  "  ds  ■  •  K: 


J:  *> 


(0)  G (6)  cos  e  de  (2.29 


This  last  integral  is  independent  of  the  radius  r,  it 
contains  only  known  functions  of  6  and  can  be  easily 
computed  either  analytically  or  numerically.  It  also 
contains  as  a  multiplicative  constant  the  stress 
intensity  factor  Kj.  This  is  fundamental  to  the 
extraction  technique.  By  choosing  the  auxiliary  functions 
<j>  and  ijj  to  be  of  the  proper  asymptotic  behavior  in  the 
neighborhood  of  the  crack  tip,  expressions  containing 
the  radius  r  disappear.  By  using  the  known  asymptotic 


expansions  of  the  displacements  in  the  neighborhood  of 


the  crack  tip,  the  stress  intensity  factor  K  appears 
as  a  multiplicative  constant.  All  the  terms  in  the 
integrals  of  and  1^  behave  exactly  in  the  same  way. 
Each  one  will  give  an  integral  of  a  function  of  9 
only,  containing  as  a  multiplicative  constant. 

Let  us  now  examine  the  integrals  of  I5  and  Ig  on 
the  arc  r  .  A  typical  term  will  be: 

S 

Substituting: 

I 

|s  .  Kj  r"1/2  Gj^e),  nx  »  -cose,  *  -  r'1/2  t(9), 

ds  =  rd6 


we  obtain: 

7 r 

G.(0)  *(9)  cos  0  d0 .  (2.30) 

-IT  1 

This  is  again  independent  of  the  radius  r  and  the  same 
arguments  apply.  It  is  then  seen  that  contour  integration 
around  the  circular  arc  will  yield  the  stress  intensity 
factor  K_  multiplied  by  a  constant. 


)  H  nx  * ds  *  -Ki  { 


In  view  of  the  singularity  of  the  auxiliary  functions 
4>  and  tp  and  their  derivatives  at  the  crack  tip,  the  exis¬ 
tence  of  the  area  integral  over  the  domain  and  the  contour 
integrals  over  r  in  (2.7)  must  now  be  examined.  Since  the 

integrands  1^  and  l2  are  highly  singular  (they  are  of 
-5/2 

order  r  in  the  neighborhood  of  the  crack  tip) ,  the 

auxiliary  functions  4  and  \p  sure  chosen  to  satisfy  the 

equations  of  equilibrium  in  that  region,  so  that  1^  and 

I2  vanish  there.  As  far  as  the  contour  integration  is 

concerned,  only  the  upper  and  lower  faces  of  the  crack 

are  of  interest.  We  normally  assume  no  tractions  there 

and  it  can  be  seen  from  (2.14)  and  (2.15)  that  Ig  and  Ig 

vanish  there.  In  the  case  where  the  crack  faces  have 

applied  tractions  on  them,  these  integrands  are  of  order 
-1/2 

r  .  The  integrals  of  Ig  and  Ig  exist,  but  care  must 
be  exercised  in  their  numerical  evaluation. 

The  integral  of  I 3  and  1^  on  the  crack  faces  requires 
further  consideration.  By  comparing  with  (2.3)  and  (2.4) 
it  can  be  seen  that  (2.10)  and  (2.11)  can  be  rewritten  as 


X(<P,  >P)u 

(2.31 

Y(<J>,  Wv 

(2.32 

where  X(4>,<P)  and  Y(<p,ty)  are  the  tractions  corresponding  to 
displacements  $  and  ip.  We  have  already  seen  that  the 


auxiliary  functions  <J>  and  \p  satisfy  the  equilibrium 

equations  in  the  neighborhood  of  the  crack  tip.  By 

considering  their  form  in  (2.25)  and  (2.26)  and 

recalling  that  the  displacement  expansions  in  the 

neighborhood  the  crack  tip  contain  terms  of  the  order 
-1/2 

r  '  we  conclude  that  the  auxiliary  functions  <P  and 
4;  can  be  chosen  to  be  the  eigenfunctions  corresponding 
to  the  eigenvalue  -1/2,  therefore  they  satisfy  the 
traction  free  boundary  conditions  on  the  crack  surfaces. 
In  other  words  X(<t>,  40  and  Y(<J>,  40  and  consequently 
and  I4  vanish  there. 

Let  us  now  summarize  the  conditions  that  the 
auxiliary  functions  <p  and  \p  must  fulfill  and  then 
proceed  to  construct  them; 

i)  4>  and  \p  must  satisfy  the  equations  of  equili¬ 
brium  in  the  neighborhood  of  the  crack  tip; 

iil  they  must  satisfy  the  traction  free  boundary 
conditions  on  the  faces  of  the  crack  near  the  crack  tip; 

iii)  they  must  have  a  singularity  of  the  order 

-1/2 

r  '  in  the  neighborhood  of  the  crack  tip; 

iv)  they  must  vanish  on  the  part  of  the  boundary 
where  displacements  are  prescribed. 

Conditions  i) ,  ii)  and  iii)  together  mean  that 
<t>  and  4;  are  eigenfunctions  for  the  problem  of  an 


infinite  body  containing  a  crack,  and  correspond  to 
the  eigenvalue  equal  to  -1/2. 

We  now  observe  that  in  the  absence  of  surface 
tractions  both  the  equilibrium  equations  (2.1),  (2.2) 
and  the  stress  boundary  conditions  (2.3)  and  (2.4)  are 
homogeneous  and  contains  only  material  constants  and 
the  derivatives  of  the  displacement  vector  components. 

By  differentiating  them  with  respect  to  x  it  is  seen 
that  ~  and  also  satisfy  the  same  equations.  There¬ 
fore  they  are  possible  candidates  for  <J>  and  satisfying 
requirements  i)  and  ii)  above.  By  comparison  of  (2.20), 
(2.21)  with  (2.25),  (2.26)  it  is  seen  that  ~  and 
also  satisfy  requirement  iii)  above.  If  we  new  consider 
the  case  of  a  body  where  only  tractions  are  specified 
on  the  contour,  then  requirement  iv)  does  not  have  to 
be  satisfied  and  we  have  arrived  at  an  explicit  form 
for  <j>  and  ip  satisfying  all  the  requirements.  By 
differentiating  (2.16)  and  (2.17)  with  respect  to  x 
and  setting  the  constant  multipliers  equal  to  one,  we 
obtain  for  mode  I: 

0  -  r"1/2  cos  |  (2-2v-5cos2  |  +4cos4  |)  (2. 


4)  *  r”^2  sin  j  (-2+2v-3cos2  j  +4cos4  f-  ) 


(2. 


and  for  mode  II: 


4>  -  r“1/2  sin  |  (-2+2v+3cos2  |  -4cos4  |  )  (2.35! 

ip  «  r  ^'//2  cos  j  (2v-5cos2  +4cos4  j  ).  (2.36! 


The  contour  integrals  on  the  circular  arc  that  appear 
in  (2.7)  can  now  be  computed  explicitly.  Omitting  inter¬ 
mediate  results,  we  obtain: 


-S 


(I3+I4)ds  + 


1 


(I5+Ig)ds 


4  (1-v)  ir  (2tt)  1/2K 


I* 

(2.37! 


Taking  into  account  (2.14)  and  (2.15)  we  obtain  from 
(2.7)  the  extraction  formula: 


4  (1-v)  tt  ( 2tt ) 


1/2 


(I^+I^dA  + 


J  (X3-fI4)ds  - 

r 


(2.38: 


In  the  case  where  only  tractions  are  specified  on  the 
boundary  of  the  body  this  simplifies  to: 


i 
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4  (l-v)  TT  (2 it) 


1/2 


K, 


■! 


(I3+I4)ds 


^  Y4>ds 
T  (2.39) 


In  the  case  of  mode  II,  the  integrals  on  the  circular  arc 
yield  the  same  constant  and  we  obtain: 


4  (l-v)  it  (2tt)  1/2  KIJ: 


(I1+I2)dA  + 


! 


(I3+I4)ds  - 


(2.40) 


and  in  the  case  of  only  applied  tractions  specified: 

4(1-v)tt(2tt)1/2Kt_  *  ^  (I3+I4)ds  -  ^  X4»ds  -  ^  Y^ds. 


r  (2.41) 


In  the  case  where  the  domain  contains  more  than  one 
geometric  singularity,  we  simply  remove  a  small  disk  of 
radius  r  ■  e  from  every  singular  point.  This  is  done 
in  order  to  be  able  to  perform  integration  by  parts. 

The  asymptotic  expansions  (2.16)  and  (2.17)  or  (2.18) 
and  (2.19)  are  considered  only  in  the  neighborhood  of 
the  point  where  extration  of  the  stress  intensity  factor 
is  desired.  The  auxiliary  functions  0  and  \p  are  also 
referenced  to  that  point.  The  other  singular  points 
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are  thus  treated  as  the  rest  of  the  boundary  and  expres¬ 
sions  (2.38),  (2.39),  (2.40)  and  (2.41)  are  still  valid. 

In  the  formulation  described  above  only  one  mode 
was  considered,  that  is,  the  vicinity  of  the  crack  tip 
was  assumed  to  experience  either  mode  I  or  mode  II 
deformation,  and  the  corresponding  stress  intensity 
factor  was  extracted.  Of  course,  in  practical  situations 
both  modes  will  be  present  and  the  two  intensity  factors 
must  be  computed  separately.  One  of  the  major  advantages 
of  the  generalized  influence  function  method  is  the 
ability  to  obtain  the  intensity  factors  for  the  two  modes 
separately.  It  will  now  be  shown  that  the  extraction 
formulae  (2.38),  (2.39),  (2.40)  and  (2.41)  are  still 
valid  when  both  modes  sure  present  simultaneously.  Let  us 
first  formally  introduce  the  following  operators  that  give 
the  tractions  in  the  x  and  y  directions  corresponding  to 
displacements  u  and  v  through  the  strain-displacement, 
stress-strain  smd  Cauchy  boundary  relations: 


X  (u,v) 


X( 


9u  3v 
3x  Ty 


+  G  ( 


3u 

3y 


+ 


+  G  ( 


.  3v  . 
+  n  ) 
3x  y 


(2.42) 


.1 

A 


i 
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Y(u,v) 


3  v 
3y 


)  n  +  G(  n  +  n  ) 
y  3x  x  3y  y' 


+ 


,  r,  9u  _  ,  9v  . 

+  G(  t-—  n  +  ir—  n  ) 
3y  x  3y  y' 


(2.43) 


Equation  (2.7)  can  now  be  written  as: 


(Ij+I^dA  - 


[X(4),^)u  +  Y(<|>,i|;)v]ds  + 


+ 


[X(u,v)$  +  Y(u,v)i|;]ds 


0. 


(2.44) 


The  only  difference  now  is  the  asymptotic  expansion  for 
the  displacements  which  is  valid  along  the  circular  arc 
Let  us  use  the  notation 

u  «  Uj  +  uXI  (2.45) 


v  ■  »i +  vn  (J-46) 

where  subscripts  I  and  II  refer  to  the  corresponding 
modes  and  Uj,  Vj,  Uj^  and  v^  are  given  by  (2.16) , 

(2.17),  (2.18)  and  (2.19)  respectively.  They  are  listed 
below  together  with  their  derivatives  as  well  as  the 
auxiliary  functions  and  their  derivatives. 
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*  a^  r^2  cos  y  (2-2v-cos2  ^  )  (2.47) 

Vj  *  dj  r^2  sin  ^  (2-2v-cos2  ^  )  (2.48) 

*  j  aj  r”^2  cos  .j  (2-2v-5cos2  +4cos4  y  )  (2.49) 

»  i  a_  r"1/2  sin  |  (2-2v-3cos2  |  +4cos4  |  )  (2.50) 

3y  2  I  2  2  2 

»  i  aT  r“1/2  sin  |  (-2+2v-3cos2  i  +4cos4  |  ) (2.51) 

dX  2  I  2  2  2 

>  j  ax  r-1^2  cos  |r  (-2v+5cos2  j  -4cos4  )  (2.52) 

*  r  ^2  cos  ^  (2-2v-5cos2  y  +4cos4  )  (2.53) 

U,  -  r”1/2  sin  |  (-2+2v-3cos2  |  +4cos4  |  )  (2.54) 

»  y  r-3^2  cos  [  6-6v+ (-43+8v)  cos2  j  + 

+  84cos4  j  -48cos^  j  ]  (2.55) 
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=  i  r"3//2  sin  *■  [2-2v+ (-23+8v)  cos2  j  + 

3y  2  2  * 

+  60cos4  y  -48cos^  j  ] 

>  i  r”3/2  sin  |  [-2+2v- (7+8v) cos2  %  + 

3x  2  2  2 

+  60cos4  j  -48cos^  j  ] 

-  I  r-3'2  cos  |  [-6v+ (35+8v) cos2  f  - 

-  84cos4  2  +  48c°s6  |-] 
Ujj  ■  ajj  r3-^2  sin  y  (2-2v+cos2  ) 

vXI  «  r1//2  cos  2  (2v-cos2  |  ) 

-  2  an  r”1/2  sin  |  [-2+2v+3cos2  j  -4cos4  ■ 


3uII  -  1  a  -”1/2 
”17”  2  aH  r 


j  [4-2v-5cos2  j  +4cos4  j 


(2.56) 

(2.57) 

(2.58) 

(2.59) 

(2.60) 

:  ] 

(2.61) 

] 

(2.62) 


cos 
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*  =  r“1//2  sin  4  [-2+2v+3cos2  f-  -4cos4  j  ]  (2.65) 

II  2  z  *• 

4;  ^  *  r”1//2  cos  j  [2v-5cos2  j  +4cos4  j  ]  (2.66) 

=  I  r"3/2  sin  |  [-2+2v+  (23-8v)  cos2  |  - 

3x  2  2  ^ 

-  60cos4  +48cos^  ]  (2.67) 

«  I  r"3/2  cos  |  [l2-6v+ (-51+8v) cos2  |  + 

+  84cos4  j  -48cos^  j  ]  (2.68) 

jyll  =  i.  r”3//2  cos  i  [6v- (35+8v)  cos2  f-  + 

dx  2  2  * 

+  84  cos4  j  -48cos®  j  ]  (2.69) 


' vn  i  -3/2  .  e  r,  .  2  e 

j  x  sxn  J  L  2v-  (15+8  v)  cos  y  + 


+  60cos4  |  -48cos6  |  ]  (2.70) 


al  = 


G/Tn 


G/I7 


(2.71) 


The  contour  integral  on  the  circular  arc  r  in 
(2.44)  can  now  be  written  as: 


-  ^  [XU,^)^  +  Y(0,^)vI]ds  - 

r‘  C. 

-  j[x(0,^)ui;r  +  Y(0,ip)vII]c 

re 

+  ^  [X(uI>vI)«J)  +  Y(uI>  Vj)  i|;]ds  + 


Y)  vn^ds 


+  \  [X(uII#  vIX)0  +  Y(uII,vIJ)^]ds. 


(2.72) 


In  order  to  extract  the  mode  I  stress  intensity  factor  we 
substitute  $  *  and  It  can  now  be  seen  that  <p , 

3$  3ti;  3uII  3vII  . 

Jt  ■  17  -  VH>  Tj-  -  -53T  *nd  nx  ■  "cose  are  sy™"etrl= 
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with  respect  to  the  x  axis  whereas  P ,  4^-  ,  1^-  ,  u__, 

o y  o x  xx 


3y  '  3x 


and  =  -sin0  are  antisymmetric  with 


respect  to  the  x  axis.  By  virtue  of  (3.42)  and  (3.43) 
we  conclude  that  XU,  ip),  X(u^,  v^)  and  Y(u^,  v  )  are 
symmetric  whereas  Y($,  ip)  ,  Y(u^,  Vj)  and  XU^,  v  ) 
are  antisymmetric  with  respect  to  the  x  axis.  Since  the 
expressions  in  the  radius  r  disappear  when  integrating  on 
the  circular  arc  r  and  the  domain  of  integration  (-it,  it) 
is  symmetric  with  respect  to  the  x  axis  it  follows  that 
the  integrals  of  XU,  ijOu^,  YU,  ijOv^,  X(uXI,  v  )  <|) , 
YCUjj,  vanish,  these  products  being  antisymmetric 

with  respect  to  the  x  axis. 

The  remaining  integrals  are: 


[XU,  ifOuj  +  Y($,  ilOVjlds  + 


^  [XtUj,  vx 


)  <t>  +  Y  (Uj,  Vj)  ^]ds 


(2.73) 


and  these  are  the  same  as  those  considered  in  the  deriva¬ 


tion  of  the  extraction  formula  for  Kj. 


In  order  to  extract  we  substitute  4>  =  $  , 


^  »  P. 


'  fy  '  fx  are  the  sYrnmetric  terms 


In  this  case  p 


-51- 


3.  EIGENVALUES  AND  EIGENFUNCTIONS  FOR  AN  INFINITE  NOTCH 

OF  ARBITRARY  SOLID  ANGLE 


Let  us  consider  a  two  dimensional  elastic  body  of 
infinite  dimensions  with  a  notch,  the  solid  angle  of 
which  is  equal  to  2a  as  shown  in  figure  (3.1).  We  choose 
the  vertex  of  the  angle  to  be  the  center  of  the  coordinate 
system.  The  bisector  line  of  the  solid  angle  is  selected 
as  the  x  axis  of  a  Cartesian  system  (x,  y)  and  the 
reference  (0=0)  axis  of  the  corresponding  polar  system 
(r,  0)  as  shown  in  the  figure.  The  faces  of  the  angle 
are  assumed  to  be  free  of  tractions.  We  shall  investigate 
the  state  of  stress  and  deformation  in  the  neighborhood  of 
the  apex.  In  the  interest  of  completeness,  all  required 
equations  are  derived  from  first  principles  [15,  16]. 

Let  us  write  the  equations  of  equilibrium  in  the  case 
where  body  forces  are  absent  in  the  form: 


(3.1) 

(3.2) 


In  the  usual  way,  let  the  stresses  be  defined  in  terms  of 
the  Airy  stress  function  x  as: 


x(ul'vl) 


(2.127) 


=  r“1/2  x1(0) 

Y(u1,v1)  -  r“1/2  Y1(0).  (2.128) 

this  integral  is  equal  to: 

'rl/2  J-TT  IX**(9)G1(9)  +  IMHl(9)ld9  + 

+  r1/2  [  U1(e)»(e)  +  Y1(e)f  (9)]de 

■  r“1/2  C5  +  r“1/2  C6  -  r“1/2  C?  (2.129) 


The  necessary  and  sufficient  condition  for  this  to  have 
the  same  value  for  all  values  of  r  is  that  C7  *  0,  that 
is  *  0. 


[X($/^)u.  +  Y  (<j)  ,\p)  v.,  ]ds  + 


r  +r  * 

1 1  3 


^  tX(u1,v1)<|>  +  Y(u1,v1 


)^]ds  =  0  (2.125) 


r  +r  * 

1  3 


therefore : 


[X  (<t> , ip) \a1  +  Y  (4> , ty)  vl]  ds  + 


)  tj;]ds 


+  ^  [X(u1,v1)<|>  + 

r3 

-  ^  [X(0,tj;)u1  +  Y(4>,^)  v1]ds  + 


^  [X(u1,v1)4>  +  Y(u1,v1)4']< 


(2.126) 


which  means  that  the  contour  integral  has  the  same 
value  on  all  circular  arcs  r£.  (In  fact  is  path 
independent  but  it  is  convenient  to  consider  only 
circular  arcs  here).  By  using  (2.112),  (2.113),  (2.77) 
(.2.80),  (2.99),  (2.100)  and  noting  that 
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u 


[Lx(<#>/I^)u1  +  Ly  (<&,  VO  v^]  dA  - 


-  ^  [X^/tjOi^  +  Y  (V>,V0  v^ds 

r 

+  ^  [Xtu^v^#  +  Y(\i1,vl 


)\|>]ds  *  0. 


(2.123) 


The  auxiliary  functions  4>  and  V>  satisfy  the  equilibrium 
equations  in  the  locality  of  the  crack  tip,  therefore 
the  area  integral  vanishes.  Both  the  auxiliary  functions 
and  the  eigenfunctions  satisfy  the  traction  free  conditions 
on  the  crack  surfaces,  therefore  the  contour  integrals 
vanish  on  the  segments  r2  an d  of  the  boundary.  Equation 

(2.123)  then  reduces  to: 


^  [X(4>,ii»)u1 


+  Y  (t^VOv^ds  + 


r  +r 

13 


^  IXCu1,v1)<()  +  Y(u1,v1)^]ds  »  0. 


(2.124) 


r  +r 

13 


We  now  choose  a  larger  annulus  bounded  by  the  circular 
arcs  r  and  r3*  and  apply  the  same  arguments.  Then: 


^  [X(u3,v3)4>  +  Y(u3,v3)i|<]ds 


.1/2 


I: 


.1/2 


[x3(e)*(0)  +  Y3(e)¥3(6) ]de  -  r  c4. 

(2.121) 


It  can  be  easily  seen  that  in  the  limit  as  t  +  0  both 
expressions  vanish. 

The  integral  L1  appears  to  be  more  troublesome  as 
it  will  be  of  the  form  r~1//2  C  and  this  is  unbounded  as 
r  ■*»  0.  We  shall  now  prove  that  this  integral  is  zero. 

Let  us  consider  a  two  dimensional  body  in  the  form 
of  am  annulus  with  a  slit,  centered  at  the  crack  tip  as 
shown  in  figure  (2.4).  The  boundary  of  this  body  consists 
of  the  two  circular  arcs  amd  r3  and  the  straight 
segments,  r2  amd  on  the  crack  surfaces.  Since  the 
eigenfunctions  satisfy  the  equilibrium  equations  we  cam 
write 


Lx(ul'vl)<j)  +  Ly (ul,vl) ^  “  0 


(2.122) 


integrating  over  the  domain  and  applying  the  Gauss 
theorem  twice  we  obtain: 
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S 


[X(u2,v2)4>  +  Y(u2,v2)t|>]ds 


I: 


[X2(0)*(0)  +  Y2(0)  4»  (© )  ]  d6  -  C2  (2.117) 


which  is  again  independent  of  r  and  cam  be  readily 
computed. 

In  order  to  evaluate  the  integral  L3  on  we  use 
(2.79) ,  (2.82)  and  observe  that: 


X (u3 ,  v3)  »  r1/2  X3(0) 


Y(u3,  v3)  -  r1/2  Y3(9) 


(2.118) 


(2.119) 


therefore : 


S 


IX  (4>,  U;)  u3  +  Y(<|),^)v3]ds 


.1/2 


1: 


1/2/ 


I%^e)G3(0)  +  Y^(9)H3(e)ld6  -  r  C3 


(2.120) 


r  *\ 


and: 
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Let  us  now  examine  each  of  these  expressions  separately. 
The  integral  L2  will  provide  the  required  K2*  In  view 
of  (2.42),  (2.43),  (2.78),  (2.81),  (2.99),  (2.100)  we 
can  write 


X(*,!|»)  * 

r"2  XM<0> 

(2.112) 

?(♦,¥)  - 

r‘2  yw(e) 

(2.113) 

X(u2»v2) 

-  r'1  X2(9) 

(2.114) 

Y(u2,v2) 

-  r"1  Y2(0) 

(2.115) 

This  last  integral  is  independent  of  r  and  can  be 
computed  either  analytically  or  numerically.  Similarly: 


where 


[X(4>,^)u  +  Y(4>»^)  v]ds  + 

•5' 

r 

X(u,v)<J>  +  Y(u,v)^]ds  *  +  Lj  +  L3 

(2.108) 

e 

• 

-h  \ 

rc 

[X(4>,ii>)u.  +  Y  ( (ji ,  ^ )  v  1  ]  ds  + 

1  x  1 

# 

+  K1 

r 

m 

\  [X (u^, v^)  <|>  +  Y(u1,v1)^]ds 

e 

(2.109) 

-k2  ) 

re 

[X  (<D ,  ip)  \x2  +  Y(<j>,4>)v2]ds  + 

0 

+ 

* 

to 

*  -1 

\  [X(u2,v2M>  +  YCu2,v2)^]ds 

e 

(2.110) 

’K3  ) 

F« 

[X(<)),^)u3  +  Y  ($ ,  41)  v^]  ds  + 

+  K, 

\  [X(u,,v.J  <f>  +  Y(u,,v.)^]ds 

(2.111) 
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Let  us  further  introduce  the  notation  L  (u,  v)  and 
Ly(u,  v)  for  the  differential  operators  of  equilibrium 
in  the  x  and  y  directions: 

Lx(u.v)  -  (A+G)  &  (£+£>  +  <2-1051 

2  2 

Ly(.u»v)  *  (X+G)  (  -|^  +  -|^  )  +  G(  j-  +  — j  )  .  (2.106) 


Equation  (2.44)  can  now  be  written: 


U 


[L  ($,i|0u  +  L  (4>,v)v]dA 
x  y 


-s 


[X(4>,4>)u  +  Y(<J>,\J;)v]ds 


(2.107) 


[X(u,v)<|>  +  Y(u,v)\f>]ds  *  0. 


At  this  point  we  are  interested  only  in  the  contour 
integrals  around  the  circular  arc  r  .  The  other  integrals 
will  not  be  affected  as  the  asymptotic  expansions  are  not 
valid  there.  The  integrals  on  r_  can  be  written  as: 
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v31(r,0)  -  r1/2  H31(0) 


(2.97) 


v32(r,0)  -  r1/2  H32(0) 


(2.98) 


The  intensity  factors  and  the  eigenfunctions  u^,  vi 
may  correspond  to  either  mode  I  or  mode  II.  The  intensity 
factor  K2  corresponding  to  the  second  term  in  the  expansion 
will  now  be  extracted.  To  this  end  the  auxiliary  functions 
are  chosen  to  be  of  the  form: 


<j>  -  r"1  4(0) 


ip  *  r_1  Y(d). 


Their  derivatives  can  be  written  in  the  form: 


(2.99) 


(2.100) 


r~2  41(0) 


r~2  $2  (0 ) 


r"2  4'1  ( 0 ) 


(2.101) 


(2.102) 


(2.103) 


r“2  ^2  (0)  , 


(2.104) 
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u1(r,9)  *  r1//2  G^t©)  (2.77) 

u2(r,0)  -  r  G2(6)  (2.78) 

u3(r,0)  ■  r3//2  G2(9)  (2.79) 

v^r,©)  «  r1/2  1^(0)  (2.80) 

v2(r ,0)  »  r  H2(0)  (2.81) 

v3(r,0)  -  r3/2  H3(9).  (2.82) 


Their  derivatives  can  be  written  in  the  form: 

v* 
*  • 

|H  -  Ki  uu(r#e)  +  K,  u21(.r,0)  +  K3  u31(r,0)  +  ... 

(2.83)  T 


-  Ki  u12(r,0)  +  K2  u22(r,0)  +  K3  u32(r,0)  +  ... 

(2.84) 


li  “  K1  vu(r'0)  +  K2  v21(r'9)  +  K3  v31(r'e)  + 

(2.85) 


and  $, 


,  ~  are  the  antisymmetric  terms.  It  follows 

that  X{$,  4>)  is  now  antisymmetric  and  Y  (<p,  40  is  symme¬ 
tric.  The  integrands  that  now  cancel  are:  X(4>,  40  u^, 

Y  C4> ,  'P)vJ,  X(ux,  vx)<p  and  YCUj,  VjH.  The  remaining 
integrals  are: 


<p)uXI  +  Y  (  <t> , 


i^)Vn]ds  + 


S 


X(u 


II 


VII}* 


+  Y  (u 


II' 


vZI)  40  ds 


(2.74) 


which  are  precisely  those  considered  in  deriving  the 
formula  for  K^. 

Another  advantage  of  the  extraction  technique  is 
the  ability  to  extract  the  intensity  factors  corresponding 
to  higher  order  terms  in  the  asymptotic  expansions.  Let 
us  write  the  asymptotic  expansions  for  the  displacements 
in  the  form: 


u  *  K!  ux(r,  9)  +  K2  u2 (r,  e)  +  K3  u3(r,  0)  +  ...  (2.75) 


v 


K1  vx(r. 


9)  +  K2  v2 (r,  9)  +  K3  v3 (r,  9)  +  ...  (2.76) 


where : 


It  can  be  easily  verified  that  the  equilibrium  equations 
are  identically  satisfied.  The  stress-strain  relations 
in  the  case  of  plane  strain  can  be  written  in  terms  of 
the  displacements  as: 


a 


x 


X( 


3u  3v 
3x  3y 


“IS 


(3. 


a 


y 


A( 


3u  3v 
3x  3y 


2 r 

2G  3y 


(3. 


G  ( 


3u  3v 
"oy  3x 


(3. 


Adding  (3.6)  and  (3.7)  yields: 


*x  +  oy-  2U  +0)  <  (3‘ 

Combining  (3.6)  and  (3.9)  and  putting  2 ( y 
obtain: 


=  v  we 


Analogously: 


2G  tt-  =  ct  -  v  (a  +  n  ) 
9y  y  x  y' 


(3.1]) 


In  the  case  of  plane  stress  we  can  write; 


a 


x 


A( 


3u  ,  9v  ,  9w 

*  TI 


'  (3.12) 


a 


y 


x  ( 


9u  .  9v  9w 
9x  9y  9z 


+  2G  1^- 
3y 


(3.13) 


0 


X( 


9u  9v  .  9w  . 
9x  3y  9z  ; 


+ 


(3.14) 


from  which  by  summation: 


ax  +  ny  *  (3A  +  2G)  (  +  f|f  ) 

and  by  substituting  back  for  +  ||r  and  putting 

3X+2G  "  T+v  we  obtain: 


2G 


9u 

7x 


+  cr  ) 

y 


(3.16) 


2G 


9v 

9y 


(ax  + 


V 


(3.17) 


We  define  a  as: 


for  plane  strain 


for  plane  stress 


(3.18) 


and  combine  (3.10),  (3.11),  (3.16)  and  (3.17)  to  obtain: 


2G  3*  =  °x  -  o(0x  +  V 


(3.19) 


2G  37  "  ay  “  a(ax  +  ay)‘ 


(3.20) 


We  now  differentiate  (3.19)  twice  with  respect  to  y, 
differentiate  (3.20)  twice  with  respect  to  x,  add  them  and 
use  (3.8)  to  obtain: 


2 

2  L'fo  =  +  ^2-  "  °v2  (°v  +  O  • 

3x3y  Ty2  3x2  X  y 


(3.21) 


Finally,  substituting  (3.3),  (3.4),  (3.5)  into  (3.21)  and 
dividing  by  (1  -  a)  we  obtain: 


V  X  -  0. 


(3.22; 


We  now  write: 


This  equation  defines  the  function  in  terms  of  the 
stress  function  x  up  to  the  functions  of  integration. 
Equations  (3.19)  and  (3.20)  can  be  written  in  the  form: 


2C  §£  -  -n  +  (1  -  a) (ox  +  ny)  (3.24) 

2G  “  ‘aX  +  t1  -  <aX  +  ay>  0.25) 

and  by  use  of  (3.23)  they  become: 

2G  !*  “  ’  3^  +  (1  “  a)  (3‘26) 

2G  ^  +  (1  *  a)  -(3*27) 

These  can  now  be  integrated  to  give: 

2G  “  =  -  +  (l-o)  ||  (3.28) 

2G  v  -  -  |i  +  (1  -  0)  ||  (3.29) 


where  u  and  v  are  known  up  to  possible  rigid  body  motion. 

These  displacements  u  and  v  were  derived  on  the 
basis  of  equations  (3.6)  and  (3.7).  The  third  equation 
of  equilibrium  (3.8)  also  has  to  be  satisfied  which 
imposes  an  one  more  condition  that  the  function  ip  has  to 
satisfy.  From  (3.8)  and  (3.5): 


Differentiating  (3.28)  and  (3.29)  and  adding: 


2  T 


xy 


(1  -  a)  72i^ 


and  in  view  of  (3.30): 
V2^  *  0. 


(3.31) 


(3.32) 


The  question  now  arises  whether  it  is  always  possible  to 
find  a  function  \ p  which  is  harmonic  and  related  to  the 
stress  function  by  (3.23).  By  taking  Laplacians  of  both 
sides  of  (3.23)  we  obtain: 


V4x 


2  2 
V  (V  X) 


3x3y 


(V2ip) 


0. 


(3.33) 


We  now  separate  the  terms  that  are  independent  of  x  and 
the  terms  that  are  independent  of  y  by  writing  <p  as : 


^  -  f^x)  +  f2(y)  +  w(x,  y)  (3.34) 

from  which: 

2  2 

V2^  «  (  +  -L_.  ,  f"(x)  +  f"(y)  +  v2w(x,  y)  (3.35) 


and  by  (3.33): 


"a*> 


— - —  (V2w)  *  0 
3x3y  '  w'  u 


(3.36) 


From  which: 


VZ  C  jy  y2  w<u'  y)3  -  0 


consequently,  the  expression  in  the  bracket  is 
independent  of  x: 


jy  V2  w (x ,  y)  »  gx(y) 


Therefore,  V  w(x,  y)  is  itself  independent  of  x: 


V2  w(x,  y)  *  V2  w (y )  -  g,(y). 


Interpreting  now  the  mixed  derivative  in  the  reverse 
order : 


&  i-h  y>3  -  0 


(3.37) 


(3.38) 


(3.39) 


(3.40) 


Therefore,  the  expression  in  the  bracket  in  independent 
of  y: 


V2  w(x,  y)  *  g3  Cx) 


(3.41) 


from  which  V  w(x,  y)  is  itself  independent  of  y: 
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V2  w(x,  y)  =  V2  (x)  =  g4(x)  (3.42) 

From  (3.39),  (3.42)  and  noting  that  constant  terms  have 
been  incorporated  in  f^(x)  and  fjCy)  it  follows  that: 

V2  w(x,  y)  =  0  (3.43) 


or 

V2  OU,  y)  -  fx(x)  -  f  j  (y) ]  “  0  (3.44) 

which  means  that  the  function  can  always  be  adjusted  by 
means  of  functions  of  x  only  and  functions  of  y  only  to 
become  harmonic.  These  functions  are  the  integration 
functions  needed  to  obtain  from  x •  There  are  no  other 
requirements  that  ij;  has  to  satisfy,  therefore,  it  is 
always  possible  to  find  such  a  function. 

Stress  functions  in  polar  coordinates 

Expressions  for  the  stress  components  in  polar 
coordinates  in  terms  of  the  stress  function  x  can  be 
derived  as  follows: 

Noting  (3.3)  and  using  the  stress  transformation  law  we 
obtain: 
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2  2 
o  «  a  cos  0  -  2T  acos0sin0  +  aQsin  6 
x  r  re  e 


(3.45) 


3  v 

transforming  — into  polar  coordinates: 

3y 


.  -L. 

cos  6  ! 

*  T 

r  i 

3 

,  cos 

3r 

r 

1 

if*  +  : 

T 

r 

2  ■ 
3  0^ 

(3.46) 


+  S2S-  |i  ) 


3,13: 


+  sin2e 

3r 


(3.47) 


comparing  the  coefficients  of  like  trigonometric  terms 
we  obtain: 


(3.48) 


re  Jr  v  r  30 
->2. 


r  fe  }  "  "  F  Tefe  + 


-M 

9  3r 


(3.49) 


(3.50) 


and  it  can  be  easily  seen  that: 


°r  *  a0  *  7  x 


(3.51) 
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The  stress-strain  relations  of  two-dimensional 
elasticity  in  polar  coordinates  can  be  written  concisely 


:r  *  TS  [*°0  +  (1  *  (ar  +  a0)] 


(3.52) 


^  t-nr  +  (1  -  <0  ("r  +  "e)3 


(3.53) 


with  a  given  by  (3.18). 

We  now  introduce  the  function  ^  which  is  defined  (up  to 
the  functions  of  integration)  as: 


7  X  -  ^  (r  -5F  > 


(3.54) 


Taking  into  account  (3.50)  and  (3.52)  the  expression  for 


er  can  be  written  as: 


IT  “  1g  C’  +  (1  “  a)  Tr  (r  TF  )] 


(3.55) 


integrating  with  respect  to  r: 


ur  *  2G  +  (1  “  o)  r  TF  ] 


(3.56) 


where  ^  now  incorporates  the  integration  function, 
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In  view  of  (3.54)  and  (3.51)  we  can  write: 


3ijj. 


3  u  v  i 

°r  +  °9  =  3r  (r  TT  5 


3i|> 


1  .  S2i^ 
36  3rT6 


3T  (*1  +  r 


3i|j. 


) 


3r 
(3.57) 


and  expression  (3.53)  for  eQ  can  be  written  as: 

3u„  .  .  a  .  13 2' 

r2  36' 


3 

+  (i  + 


(3.58) 


Multiplying  through  by  2Gr  and  substituting  for  ur  from 


(3.56)  gives  after  some  cancellation  of  some  terms: 

2 


3u 


2G 


e 


30 


-  ;4+a  -  o) 

r  36^ 


d  <p. 


9r36 


(3.59) 


integrating  with  respect  to  6: 


3ty. 


u9  -  ^  7  U  +  11  -  °>  r"  TF  3 


(3.60) 


where  once  again  the  function  of  integration  has  been 


incorporated  in  ip  . 


So  far  only  the  radial  and  circumferential  components 
of  stress  and  strain  have  been  considered.  The  considera¬ 
tion  of  the  stress-strain  and  strain-displacement  relations 
for  shear  stress  and  strain  will  lead  to  an  extra  condition 
that  the  function  ^  must  satisfy.  From: 


4 

J 


.‘-I 


J 


J 


l" 
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T 


re 


(3.61) 


and: 


Y 


re 


A  *  !^e  .  ^ 

r  36  3r  r 


(3.62) 


Substituting  for  ur  and  uQ  from  (3.56)  and  (3.60)  after 
some  manipulation  we  obtain: 


r  36' 


A  in  +  !!!i  )  +  1  iz.  ] 
r  9r  17^  7  36 


or 


re 


5  t-  +  T  (1  -  0)'72*1  *  -tH  ]  <3-«> 


Comparing  with  (3.56)  and  (3.61)  we  deduce  that: 


v2  n»1  =  o 


(3.64) 


The  solution  of  the  two-dimensional  elasticity  problem 
then  reduces  to  the  determination  of  a  biharmonic 
function  x  and  a  harmonic  function  4^  related  to  it 
through  (3.54). 


0 


(3.68) 


F1V(8)  +  2  (<  2+l)  F"  ( 6 )  +  (k-1)2(<  +  1)2F(9)  = 

Let 

m9 

F (0)  *  a  e 

The  characteristic  equation  becomes: 

m4  +  2  (<2+l)m2  +  (<-l)2(<+l)2  =  0 
which  has  the  roots: 

m  *  +  (k  +  1) i 

m  *  +  (<  -  1)  i 

these  roots  leading  to  solutions  of  the  form: 

sin(x+l)9,  cos(k  +  1)9,  sin(<-l)0,  cos(<-l)9 

and  the  general  solution  can  then  be  written  in  the 
form: 

X  »  r<+1[C1sin (<+l)  9+C2cos  (<  +  l)  9+ 

+  C^sin  (<-l)  0+C4cos  (<-l)  9] 


(3.69) 


(3.70) 


(3.71) 


3.72) 
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At  this  point  it  should  be  noted  that  the  above  form  for 
X  satisfies  the  biharmonic  equation  for  any  value  of  < , 
real  or  complex. 

In  order  to  determine  ^  we  investigate  solutions 
of  the  form: 


=  rm  G  ( 0 ) 


Since  is  harmonic,  from 


( 


1  3  _1_ 

1  3r  r2  362 


we  obtain: 


G"  (0)  +  m  G(0)  *  0 


(3.73) 


(3.74) 


Let 


G  = 


then: 


2  2 
n  +  m  =  0 

from  which: 


(3.75) 


(3.76) 


n  *  +  mi 


(3.77) 


G2(r,  8).  Thus  a  complex  root  of  the  eigenequation  gives 
rise  to  two  distinct  eigenfunctions.  This  can  also  be 
seen  if  we  write  the  stress  intensity  factor  in  complex 
form  and  write  the  term  corresponding  to  (3.140)  as: 

u  —  (K1+iK2)  r  [G1(r/0)+iG2(r,8)]  (3.141) 

separating  real  and  imaginary  terms : 

<1  k  . 

u  ■  r  G1(r,8)  -  K2  r  G2(r,8)  + 

+  i  r  X[K1  G2(r,6)  +  K2  G1(r,0)]  (3.142) 

Both  the  real  and  the  imaginary  parts  of  this  expression 

are  acceptable  solutions.  However,  we  observe  that  the 

imaginary  part  contains  the  same  eigenfunctions  G^(r,  8) 

G2(r,  0)  as  the  real  part.  Therefore,  to  a  complex 

eigenvalue  there  correspond  two  distinct  eigenfunctions, 

each  with  its  own  stress  intensity  factor.  They  have  a 

K1 

singularity  of  the  same  strength  r  .  If  the  complex 

number  in  (3.136)  is  an  eigenvalue,  its  complex  conjugate 

will  also  be  an  eigenvalue.  To  this  eigenvalue  there 

correspond  two  distinct  eigenfunctions  G2(r,  0),  G4(r,  8), 

<1 

both  having  a  singularity  of  the  same  strength  r  .  We 
then  conclude  that  to  a  pair  of  complex  conjugate 
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It  is  seen  that  only  the  real  part  of  the  complex  root 
determines  the  strength  of  the  singularity,  the  imaginary 
part  being  incorporated  in  the  eigenequation.  The 
expressions  in  the  bracket  in  the  eigenequations  (3.130), 
(3.121)  of  mode  I  and  (3.130),  (3.131)  of  mode  II  will 
also  be  complex.  Let  us  write  (3.120)  for  example  in 
the  form: 

u  *  r<[F1(9)  +  i  F2(0J]  (3.138) 

In  view  of  (3.148)  this  can  be  written  as: 
u  *  r  (9)  cos  (<2Jlnr) -F2  (0)  sin(<2£nr) 

+  i(F2(9)cos  (<2Jlnr)  4^(0)  sin  (<2i.nr)  )]  (3.139) 

The  terms  cos(<2inr)  and  sin(*2inr)  are  highly  oscillatory 
in  the  neighborhood  of  the  notch  tip,  the  frequency  of 
oscillations  tending  to  infinity  as  the  radius  r 
approaches  zero.  Let  us  now  write  (3.139)  in  the  form: 

<  i 

u  ■  r  [G1(r,0)  +  i  G2(r,0)]  (3.140) 

There  will  be  a  stress  intensity  factor  corresponding 
to  the  real  part  G^(r,  0)  of  the  eigenfunction  and  a 
different  one  corresponding  to  the  imaginary  part 


-78- 


2GC  (2it)  ^■//2v  =  K*  r<  [  (CB+<  +  l)  sin<  9-2KCsin0cos  (<-l)  6] 

(3.133) 

and  for  mode  II: 

2GC  (2tt)  ^^u  =  K*^  r<[  (CD-<-l)  cos<6+2<Csin0cos  (<-l)  0] 

(3.134) 

2GC(2it)1/2v  =  K*x  r<[- (CD+<  +  l)cos<0-2<Ccos9cos  (<-l)  0] 

(3.135) 


It  is  reminded  that  the  eigenvalue  <  and  the  constant  C 
are  different  for  each  mode. 

Behavior  of  the  roots  of  the  eigenequations ■  The  case 
of  complex  eigenvalues. 

The  above  derivation  is  valid  for  complex  values  of 
the  eigenvalue  <  as  well.  Let  us  write  the  complex  roots 
in  the  form: 

<  »  <1  +  i  <2  (3.136) 

With  this  notation,  the  singular  term  can  be  written  as: 

K  . 

r  *  r  cos(<2S'nr)  +  i  sin(<2inr)  (3.137) 


HkmSkt 


..  ^  ^  .*  ,•  k,  s  *.ji  v\.  '  ft  v  v  V  V.-.V  »r;  ^rM'  v.  5rj  r,  «  .  *■;  r-: 
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These  are  transformed  into  Cartesian  coordinates  and 

/ 

after  trigonometric  transformations  they  are  written 
in  the  form: 

2Gu  *  r<[ (CD-<-l) sin< 0+2KCsin6cos (<-l) 0]  (3.130) 

2Gv  *  rK[- (CD+<+1)cosk0-2kCcos0cos (<-l) 0]  (3.131) 

where  D  is  given  by  (3.128)  and  C  is  the  ratio  in  (3.100). 

The  eigenfunctions  (3.120),  (3.121)  of  mode  I  and 
(3.130),  (3.131)  of  mode  II  are  determined  within  a 
multiplicative  constant:  any  multiple  of  these  eigen¬ 
functions  will  also  be  an  eigenfunction.  This  multipli¬ 
cative  constant  is  chosen  to  agree  with  the  one  adopted 
in  engineering  literature  on  Linear  Elastic  Fracture 
Mechanics.  The  generalized  stress  intensity  factors, 
which  are  the  amplitudes  of  the  singular  terms  in  the 
displacement  and  stress  expansions  in  the  neighborhood 
of  the  reentrant  corner,  are  then  defined  to  coincide 
with  Irwin's  definition  [13]  of  and  by  writing 
for  mode  I: 

2GC  (2ir)  ^^u  *  K*  r <[  (CB-k-1)  cosk  0-2<Ccos0cos  (k-1)  0  ] 

(3.132) 


hhhii 
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k~J 

F(6)  =  sin (<+1)6  +  C  sin(«-l)8 

(3.122)  “ 

<3  (0)  =  ~  c  cos(k-I)  6 

»  ■ 

(3.123)  3 

from  which: 

F'CQ)  *  (<+l)cos  (<+l)  6  +  C (<-l) cos  (<-l)  6 

(3.124) 

G1  (0)  =  4C  sin(<+l)0 

(3.125) 

and  the  displacements  for  mode  II  can  be  written 

as :  V 

2Gur  *  rK[-  (<+l)  sin  (k+1)  0+CDsin  (<-1)0] 

► ' 

(3.126)  l 

2Gufl  =  r^C- (<+l)cos  (<  +  l)  9-CEcos  (<-l)  8] 

0 

(3.127) 

where : 

D  ■  4  ( 1—  O’)  -  <  -  1 

(3.128) 

E  *  4  (1-n)  +  <  -  1 

(3.129) 

■> 

‘  /  -•  /  /  .■  ,*  \*  *.*  /,  **  -*/•’,  •*’,  •**.  •*.  V.  *' »  * 

y 

(3.116) 


u  =  cos6  -  sin8  u 

v  »  sin6  u_  +  cosQ  ua  (3.117) 

i  U 

to  be: 

2Gu  *  r<[- (<+l)  cos< 6+CAcos6cos  (<-l)  6 

-  CBsinGsin (<-l) 9]  (3.118) 

2Gv  «  r<[  (<+l)  sinKe+CAsinQcos  (<-l)  9+CBcos0sin(<-l)  6] 

(3.119) 

Through  trigonometric  transf onnations  these  can  be 
written  in  the  form: 

2Gu  *  rK[ (CB-k-1) cos<0-2kCcos0cos (<-l) 0]  (3.120) 

2Gv  ■  r<[  (CB+<+l)  sin0-2<Csin0cos  (<-l)  8]  (3.121) 

where  B  is  given  by  (3.115)  and  C  is  the  ratio  in  (3.102). 
Mode  II 

In  a  similar  way  by  taking  »  1  and  calling  the 
ratio  *  c  in  (3.100)  the  stress  functions  can  be 


written  as: 


and: 


G  (0)  =  y  C  sin(<-l)0 


(3.109) 


from  which: 


P'  (9)  *  -  (k+1)  sin (k+1)  0  -  C  (<-l)  sin (<-l)  6 


(3.110) 


G'  (9)  a  4C  COS  (<-1)0 


(3.111) 


using  (3.82)  and  (3.83)  we  can  write  the  displacements 
for  mode  I  as: 


2G  u  ■  r<[- (k+1)  cos  (k+1)  0+CAcos  (<-l)  0] 
r 


(3.112) 


2G  uQ  a  r  [  (k+1)  sin  (k+1)  0+CBsin  (k-1)  0] 


(3.113) 


where: 


A  -  4(1-0)  -  (K+1) 


(3.114) 


B  -  4(1-0)  +  K-l 


(3.115) 


The  displacements  u  and  v  in  Cartesian  coordinates  can  now 
be  found  through  the  transformation  law: 
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ir 


I 

■ 


In  the  case  of  a  crack  (a  =  180°) ,  the  first 
positive  eigenvalue  <  *  1/2  is  of  great  significance  in 
linear  elastic  fracture  mechanics,  where  the  symmetric 
mode  is  termed  mode  I  and  the  antisymmetric  mode  is 
termed  mode  II.  In  the  sequel  this  terminology  will  be 
extended  to  the  case  of  a  notch  of  an  arbitrary  angle. 

In  the  crack  case  the  two  modes  have  the  same  eigenvalue 
but  the  eigenfunctions  corresponding  to  them  are 
different.  In  all  other  cases  (except  the  half  plane 
case  as  we  saw)  the  two  modes  have  different  eigenvalues. 
This  means  that  whenever  and  are  nonzero  then  C2 
and  C4  must  be  zero  and  vice-versa. 

It  is  important  to  note  that  in  either  mode  I  or 
mode  II  if  <  is  an  eigenvalue,  so  is  -<••  This  is  very 
important  for  the  formulation  of  the  generalized  influence 
function  method,  where  use  is  made  of  the  eigenfunctions 
corresponding  to  these  negative  eigenvalues.  The  negative 
roots  yield  unbounded  displacements  at  the  vertex  of  the 
angle,  so  they  do  not  have  physical  significance. 
Expressions  for  the  displacements 
Mode  I 

By  taking  C2  *  1  and  calling  the  ratio  C4/C2  *  C  in 
(3.102)  we  can  write: 


r  - 


F ( 6)  »  cos(<+l)9  +  Ccos(k-1)9 


(3.108) 


hold  simultaneously.  Excluding  the  meaningless  case  of 
a  *  0,  these  are  satisfied  either  when: 

a  *  90°  and  k  »  n  (half  plane)  (3.105) 

or: 

a  *  180°  and  <  «  j  (crack  case)  (3.106) 

with  n  =  integer. 

In  the  case  of  the  half  plane  *  0  and  is  undetermined. 
The  triviality  of  the  case  <  -  0  can  be  seen  by  expressing 
rr  in  terms  of  k.  Making  use  of  (3.50)  and  (3.72)  we 
obtain : 

2 

+  rK  ”1  [F"  (  0)  +  (k+1)  F  (  0)  ]  (3.107) 

r  r  302  r  dZ 

and  for  <  =  0  it  can  be  verified  that  the  expression  in 
the  bracket  becomes  zero.  From  (3.84)  and  (3.85)  it  is 
obvious  that  a g  and  xrQ  are  also  zero,  and  from  (3.82) 
and  (3.83)  we  can  see  that  the  displacements  do  not 
depend  on  the  radius  r,  therefore,  we  have  a  state  of 
rigid  body  motion.  It  is  also  of  interest  to  note  that 
for  the  antisymmetric  case  (equation  3.99),  the  value 
<  *  1  is  a  solution  for  any  notch  angle  a. 
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C3  _  sin  (<+l)  a  _  (<-4-1)  cos  (tc  +  1)  a  ,,  ,  f  n, 

Cx  ‘  "  sin  (k-1)  a  ”  ”  (<-l)  cos  (<-l)  a 

but  one  of  them,  say  C1,  can  be  assigned  arbitrarily. 
Conveniently  we  choose  »  1. 

Similarly  the  other  two  equations  (3.96)  and  (3.97) 
yield  the  eigenequation: 

sin2a  +  sin2<a  «  0.  (3.101) 


Again  for  every  root  of  this  equation  at  least  one  of  C2* 
C4  is  nonzero.  One  of  them  can  be  given  an  arbitrary 
value  but  their  ratio  is  determined  from  either  (3.96)  or 
(3.97)  as: 
r 

4  _  cos  (k+1)  a  _  (k+1)  sin  (<+l)  a  , 

Cj  ~  cos  (<-l)  a  “  ”  (<-l)  sin(ic-l)  a 

Examining  equations  (3.99)  and  (3.101)  together  it  is 
seen  that  they  are  satisfied  simultaneously  either  in 
the  trivial  case  <  *  0,  or  when: 


sin2<a  *  0 


(3.103) 


and: 


sin2a  *  0 


(3.104) 


(<+l)  C2sin  (<+l)  a  +  (<-l)  C^sin  (<-l)  a  =  0 


(3.97] 


It  is  then  seen  that  the  original  system  of  four  equations 
in  four  unknowns  is  separated  into  two  independent  systems , 
each  of  two  equations  in  two  unknowns.  It  will  be  seen 
that  (3.94)  and  (3.95)  lead  to  a  solution  which  is  anti¬ 
symmetric  about  the  x  axis  and  (3.96),  (3.97)  lead  to  a 
solution  which  is  symmetric  about  the  x  axis.  This 
separation  into  the  symmetric  and  the  antisymmetric  terms 
could  have  been  anticipated  considering  the  symmetry  ox 
the  domain  about  the  x  axis.  Each  system  is  homogeneous 
and  for  non-zero  solutions  to  exist  the  determinant  of 
each  one  should  vanish.  The  determinant  of  the  system 
(3.94),  (3.95)  is  zero  if: 


(k-1)  sin(<  +  l)  acos  (<-l)  a-  (<+l)  sin (<-l)  acos  (<+l)  a  =  0 


(3.98) 


this  simplifies  to: 


<sin2a  -  sin2<a  *  0. 


(3.99) 


For  every  root  <  of  this  eigenequation  at  least  one  of 
and  is  nonzero.  The  ratio  C^/C^  can  determined 


from  either  (3.94)  or  (3.95)  as: 
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sub  stitu  ting  for  F ( 0 )  from  (3.72)  yields: 

C^sin  (k+1)  ct+C2cos  (k+1)  a+C^sin  (k-1)  ot+C4cos  (k-1)  oc  *  0 

(3.90) 

-C^sin  (<+l)  a+C2cos  (<+l)  a-C^sin  (<-l)  a+C4cos  (<-l)  a  =  0 

(3.91) 

(k+1)  C1cos  (k+1)  a-  (<+l)  C2sin  (k+1)  a+  (k-1)  C3cos  (k-1)  a 

-  (k-1)  C4sin  (k-1)  a  *  0  (3.92) 

(k+1)  C1cos  (k+1)  a+  (k+1)  C2sin  (k+1)  a+  (k-1)  C3cos  (k-1)  a 

+  (k-1)  C4cos  (k-1)  a  *  0  (3.93) 

By  simple  additions  and  subtractions  these  are  easily 
transformed  into: 

C^sin(K+l)a  +  CjSinlK-Da  *  0  (3.94) 

(k+1) C1cos  (k+1)  a+ (k-1)  C3cos  (k-1)  a  =  0  (3.95) 

and: 

C2cos(<  +  l)a  +  C4cos(<-l)a  »  0  (3.96) 


V  .**  % 


A*  -  *  «  ' 


The  characteristic  values  of  <  for  which  the  stress 


function  in  (3.72)  and  consequently  the  displacements  in 
(3.82)  and  (3.83)  provide  a  solution  for  the  notch  pro¬ 
blem  will  be  determined  by  the  requirement  that  the 
traction-free  boundary  conditions  be  satisfied  on  the 
edges  of  the  angle.  Let  us  express  a0  and  xrQ  in 
terms  of  <.  Using  (3.49),  (3.50)  and  (3.72)  we  obtain: 


a.  =  =  K(<+l)r<_1  F  ( 6 ) 

6  3r 


_ 1  9 X .  .  _L  lx _ *-l  F  •  rn 

re  “  r  3r  36  f2  Te  “  <r  F  l0) 


(3.84) 

(3.85) 


In  order  for  tractions  to  be  zero  on  the  two  faces  of 
the  angle,  both  aQ  and  TrQ  must  vanish  there: 

nQ(e  =  +  a)  =0  (3.86) 


xrQ (0  *  +  a)  =  0 


(3.87) 


from  which: 
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which  leads  to  solutions  of  the  form: 

sinmQ,  cosmQ 

The  general  solution  for  G ( 0 )  will  then  be: 

G  =  a^  cosmQ  +  a2  sinmQ  (3.78) 

and  this  is  valid  for  complex  values  of  m  as  well. 

From  the  relationship  (3.54)  between  x  and  iK  by  equating 
the  powers  of  r  we  obtain: 

k  -  1  »  m  (3.79) 


and  by  equating  subsequently  coefficients  of  similar 
trigonometric  terms  we  obtain: 


al 


4 

K-l 


(3.80) 


a2 


4 

K-l 


(3.81) 


The  expressions  for  the  displacements  now  become: 


2Gur  *  r*  [-  (k+1) F  { 6)  +  (1-n)  G’  ( 6)  ]  (3.82) 

2Gu0  *  rK[-F'  (0)  +  (l-n)  (A-l)G(Q)  ]  . 


(3.83) 


eigenvalues  there  correspond  four  independent  eigen- 

<1 

functions  with  a  singularity  of  the  same  strength  r 
but  each  having  its  own  stress  intensity  factor.  By 
using  the  complex  form  of  trigonometric  functions 

sin  <  =  sin  cosh  <2  +  i  cos  sinh  <2  (3.143) 

cos  <  =  cos  cosh  <2  -  i  sin  sinh  <2  (3.144) 

and  substituting  in  (3.120),  (3.121),  (3.130),  (3.131) 
after  rather  lengthy  but  otherwise  straightforward 
trigonometric  and  algebraic  transformations  we  obtain 
the  following  explicit  expressions  for  the  eigenfunctions: 
mode  I : 

4G  (2tt)  1/2  n1  -  Kj15  r  1  ^e  (<1+3-4a-L1)  ^  + 

+  (i<2_L2^s1  "  2  (<1c3+<2s3)  cos6] 

-K,0 

+  e  [ (<1+3-4a-L1) c2  -  (<2»L2)s2  - 


2 (<n  C.+K,S.) COS0] 


(3.145) 


4G  (2ir) 


r  (e  [  (<  ^+3-4rr+L^) 


4G(2tt) 


4G(2tt) 


W*  „  = 

V1  KI 


-  (<2+L2)c1  -  2  (<1c3+<2s3)  sin0] 


~K  29 


+  e  [  (<1+3-4a+L1)  s2  +  (<2+L2)c2  + 


+  2 (-< 1c^+< 2s4) sine] 


} 


(3.146) 


1/2  u2  =  K22)  r  1^e  2  [-(<1+3-4a-L1) sx 


+  (k2”L2)c1  +  2  (<ls3“lc2C3)  cos9^ 


-<2e 


+  e  [  (<1+3-4n-L1)  s2  +  (<2"L2)c2  “ 


-  2 (<1s4+<2c4) cose] ^ 


(3.147) 


1/2 


v,  *  K 


<2)  r<1{e<26[«1+3-4a+L1)c1 


+  (<2+L2)s1  +  2  (<is3”k2c3)  sine] 


-<2e 


+  e  [- (<1+3-4n+L1)  c2  +  (<2+L2)s2  - 


2 (k  s „c. ) sin6 


(3.148) 


where 


c1  *  cos (<18-<2^nr) 

*  sin (K^0-K2&nr) 
c2  =  cos (<18+<2inr) 
s2  *  sin (<^0+<2inr) 
c^  =  cos (<^e-0-<2inr) 
s^  *  sin (<^0-0-<2Jlnr) 

p 

c.  =  cos  (< .  0-0+<-Jlnr) 
4  12 


(3.149) 

(3.150) 

(3.151) 

(3.152) 

(3.153) 

(3.154) 

(3.155) 

(3.156) 

(3.157) 

(3.158) 


(3.159) 


cosh2<  2acos2a+cos2K 
C1  cosh2<2a+cos2 (k^-1) a 

sinh2<2Ct  sin2<^a 
*"2  “  cosh2<  2a+cos2  (k  ^-1)  a 

and  a  is  given  by  (3.18) 
mode  II 


(3.160) 


4G(2ir)1/2  ux  *  K^}  r  1^e  2  [  (3-k  1-4ct-L1)  a1  + 

+  (<2+L2)c1  +  2  (<1c3+<2s3)  sine] 

0 

+  e  2  [  (3-«)c1-4n-L1)  s2- -  (<2+L2)c2  + 

+  2 (<1c4-<2s4) sine] ^  (3.161) 

4G(2tt)1/2  v1  =  K^}  r  1(e  2  [- (3-<1«4a+L1)  + 


+  (<2-L2)s1  +  2(-<1c3+<2s3)cos6] 

-<2e 

+  e  [- (3-<1+4a+L1) c2  -  (< 2"L2^  s2  “ 


2 (< , c .+< -s  . ) cos6 


(3.162) 


4G  { 2tt )  x/ ^  u2  =  Kjj  r  {e  [  (3-<1«4a+L1)  + 


+  (<2-L2)s1  +  2 (-<1s3+<2c3) sine] 


— <  0 

+  e  2  [- (3-<1-4a-L1)  c2  -  (<2-L2)s2  + 


+  2  (»c^s^+<2c4)  sine]  ) 


(3.163) 


0 

4G(2tt)1/2  v2  =  Kj2)  r  L{e  2  [  O-K^a+I^)  s1  + 


+  (<2-L2)c1  +  2  (<1s3-<2c3)cos6] 


-<2e 


+  e  [- (3-<1-4a+L1)  s2  +  (|C2*’L2)c2  - 


-  2 (k1s4+<2c4) cose] } 


(3.164) 


where 

and  s^  are  given  by  (3.149)  through  (3.155) 


C1**1+1)+<2C2 


_  Am  —  Sm 

C1  C2 


(3.165) 


C1<2"C2 


(3.166) 
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-cosh2<  2acos2a+cos2< 
C1  =  cosh2<2a-cos2 (<1“1) a 

sinh2<2a  sin2<^a 
^2  “  cosh2K2a-cos2  (k  ^-1)  a 


(3.167) 


(3.168) 


The  heuristic  approach  that  was  followed  in  this 
chapter  can  be  justified  in  the  framework  of  a  very 
general  theory  given  by  Kondratev  [17].  He  showed  that 
the  solution  of  an  elliptic  boundary  value  problem  in 
the  neighborhood  of  an  angle  can  be  written  in  the  form: 

u  =  I  ai  r  *  !nq  r  <J>i(6)  +  w(r,  0)  (3.169) 

i 


where 

g  *  0 ,  .  • .  $  ni^  *  1 


nu  *  the  multiplicity  of  the  eigenvalue  <± 
<j>^(0)  is  a  smooth  function 
and  w(r,  9)  is  smoother  than  the  r  1  £nq  r  terms 
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It  is  now  of  interest  to  explore  further  the  nature 
of  the  roots  of  the  eigenequations.  We  are  only 
interested  in  roots  with  a  positive  real  part.  The  roots 
with  a  negative  real  part  are  readily  obtained  from  the 
former  through  a  simple  sign  reversal,  due  to  the  symmetry 
of  the  eigenequations  with  respect  to  the  x  axis.  They 
are  disregarded  in  the  displacement  expansions  as  they 
would  imply  unbounded  displacements  at  the  tip.  Karp  and 
Karal  [18]  have  shown  that  there  exists  a  nontrivial  real 
root  which  is  always  smaller  than  the  positive  real  part 
of  any  of  the  complex  roots. 

The  smallest  positive  eigenvalue  for  mode  I  is  given 
in  table  3.1  for  various  values  of  the  half  angle  a. 

Table  3.2  displays  the  same  information  for  mode  II.  We 
can  see  that  except  for  the  crack  case  (a  *  180  degrees) , 
where  the  eigenvalues  for  the  two  modes  are  the  same,  the 
eigenvalue  for  mode  II  is  always  larger  than  the  eigen¬ 
value  for  mode  I.  When  the  half  solid  angle  a  is  equal  to 
128.7268  degrees  the  eigenvalue  for  mode  II  is  equal  to  1, 
therefore  the  stresses  corresponding  to  mode  II  are  bounded 
at  the  notch  tip.  Under  mode  I  conditions  any  reentrant 
corner  will  give  rise  to  singular  stresses  at  the  tip,  the 
critical  half  angle  being  90  degrees  (half  plane) .  The 
second  and  higher  eigenvalues  are  complex  for  angles  a 


equal  to  or  smaller  than  157  degrees  in  the  case  of  mode  I, 
and  for  angles  ct  equal  to  or  smaller  than  164  degrees  in 
the  case  of  mode  II.  In  the  case  of  a  crack  (a  =  180 
degrees)  all  the  eigenvalues  are  real.  For  smaller  angles 
there  is  a  finite  number  of  real  eigenvalues  and  an 
infinite  number  of  complex  ones.  Table  3.1  refers  to  mode 
I  and  displays  the  first  nine  real  eigenvalues  for  a 
equal  to  175  degrees  together  with  all  the  real  eigenvalues 
for  a  equal  to  170,  165,  160  and  158  degrees.  Table  3.4 
displays  the  corresponding  results  for  mode  II,  with  the 
first  nine  real  eigenvalues  for  a  equal  to  175  degrees 
and  all  the  real  eigenvalues  for  a  equal  to  170  and  165 
degrees . 
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Table  3.1 

Smallest  positive  eigenvalue 
for  various  solid  angles,  mode  I 
(roots  of  the  equation  <sin2a  +  sin2<a  =  0) 


solid  angle 

a 

first  eigenvalue,  mode  I 

360° 

180° 

0.5 

350° 

175° 

0.500052987126443 

340° 

170  ® 

0.500426375426056 

330° 

165° 

0.501453008713551 

320° 

160° 

0.503490483184783 

310° 

155° 

0.506932842286465 

300° 

150° 

0.512221361160512 

290° 

145° 

0.519854303113919 

280° 

140° 

0.530395719129773 

•  270° 

135° 

0.544483736782464 

260° 

130° 

0.562839480481682 

250° 

125° 

0.586278864957285 

240° 

120  0 

0.615731059490783 

230° 

115° 

0.652269555181627 

aar” 

220° 

110° 

P,6-73mfj 

0.697164972097201 

210° 

105° 

0.751974545407642 

200° 

100° 

0.818695851323838 

190° 

95° 

0.900043811488137 

180° 

90° 

1. 
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Table  3 . 2 

Smallest  positive  eigenvalue 
for  various  solid  angles,  mode  II 

(roots  of  the  equation  xsin2a  -  sin2<ct  =  0) 


solid  angle 

1 

a 

i 

first  eigenvalue,  mode  II 

360° 

180° 

0.5 

350° 

175° 

0.529354738341384 

340° 

170° 

0.562006549619481 

330° 

165° 

0.598191849614085 

320° 

160° 

0.638182471293363 

310° 

155° 

0.682294830307061 

300° 

150° 

0.730900741512950 

290° 

145° 

0.784440552974094 

280° 

140° 

0.843439568929300 

270° 

135° 

0.908529189846099 

260° 

130° 

0.980474925453105 

257.4536° 

128.7268° 

1. 

250° 

125° 

1.060214662528446 

240° 

120° 

1.148912751316944 

230° 

115° 

1.248039607030766 

ses; 

220 

110° 

1.359494953661662 

210° 

105° 

1.485811706900859 

200° 

100° 

1.630525086564494 

190° 

95° 

1.798932622346293 

180° 

90° 

2. 

1 

•1 

.1 

j 


*1 


Table  3.3 


Real  eigenvalues  for  various  solid  angles,  mode  I 
(Roots  of  the  equation  <sin2a  +  sin2<a  =  0) 


Solid  angle 

a 

real  eigenvalues,  mode  I 

350° 

175° 

0.500052987126443 

1.058842953176205 

1.499727767815282 

2.118822841754057 

2.497979910848996 

3.181532712089234 

3.493301550223319 

4.250184120916684 

4.482534709491662 

etc. 

340° 

170° 

0.500426375426056 

1.125406650991640 

1.497613486365886 

2.267186596933758 

2.476769998913093 

330° 

165° 

0.501453008713551 

1.202957173241424 

1.490377798463559 
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Table  3 . 3 
(continued) 


Solid  angle 

0L 

real  eigenvalues,  mode  I 

320° 

160° 

0.503490483184783 

1.302693359118874 

1.467008439164243 

316° 

158° 

0.504675297031383 

1.365731350131872 

1.436282183283593 

Table  3.4 


Real  eigenvalues  for  various  solid  angles,  mode  II 
(Roots  of  the  equation  <sin2a  -  sin2<a  =  0) 


Solid  angle  a  real  eigenvalues ,  mode  II 

0.529354738341384 

1. 

1.588609191187519 

1.999106964391621 

2.649698718055999 

350°  175° 

2.996140969135392 

3.714772677395547 

3.989019694329133 

4.789341976925038 

» 

4.972275638161229 

etc. 

0.562006549619481 

1. 

1.692250101505664 
340°  170°  1.991384797275372 

2.883886605832655 

2.920169039793589 

0.598191849614085 

330°  165°  lm 

1.838934252571961 

1.948555887250050 


2/2 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  of  STAN0AR0S-i963-A 


Figure  4.2 

Cartesian  and  polar  coordinate  systems 


-97 


K  K 

u  *  Kx  r  1  U1(0)  +  K2  r  2  U2(0)  +  ... 


(4.1) 


K  1C 

V  «  X  r  1  V^©)  +  K  r  2  v2(0)  +  ... 


(4.2) 


where  K^,  K2  are  the  generalized  stress  intensity  factors  and 
<2  are  the  eigenvalues.  The  precise  form  of  functions 
U±(e)  and  V^(0)  was  derived  in  the  previous  chapter.  We  are 
interested  in  compute  so  the  generalized  stress  intensity 
factors  K^,  K2,  etc.  The  eigenfunctions  U^(0),  V^(e)  also 
depend  on  the  elastic  constants,  the  half  solid  angle  a  and 
the  corresponding  eigenvalue  ic^. 

We  choose  as  extraction  functions  <p  and  ip  the  eigenfunc¬ 
tions  corresponding  to  the  negative  eigenvalue  of  the  term, 
whose  amplitude  is  being  extracted. 


4>  *  r  4>  (0)  + 


(4.3) 


\p  *  r  ¥(0)  +  \p. 


(4.4) 


where  $(.0)  and  f(e)  are  obtained  from  1^(0)  and  V1(e)  respec¬ 
tively  if  we  substitute  -<^  for  k^.  The  terms  ^  and  are 
smooth  functions  that  vanish  in  the  neighborhood  of  the  notch 
tip  and  arc  so  chosen  that  the  extraction  functions  satisfy 
the  following  requirements: 

i)  They  vanish  on  the  part  of  the  boundary  where  dis¬ 
placements  are  specified.  This  is  accomplished  by  the 
audition  of  the  smooth  functions  and  ip^,  which  are  called 


blending  functions.  Another  way  is  to  multiply  the  singular 
part  in  <j>  and  ip  by  a  cut-off  function,  which  takes  the  value 
of  unity  in  the  neighborhood  of  the  notch  tip  and  the  value 
of  zero  away  from  it  with  a  smooth  transition  in-between. 
Further  discussion  on  the  cut-off  and  blending  functions  can 
be  found  in  [3]. 

ii)  They  satisfy  the  equilibrium  equations  in  the 
neighborhood  of  the  notch  tip.  If  no  displacement  boundary 
conditions  are  specified  on  any  part  of  the  boundary,  then 
they  satisfy  equilibrium  everywhere  in  the  domain. 

iii)  They  satisfy  the  traction-free  boundary  conditions 
on  the  faces  of  the  notch,  that  is  tractions  corresponding  to 
displacements  <j>  and  \p  vanish  for  0«a  and  9*-a. 

-<1 

iv)  They  have  a  singularity  of  the  order  r 

The  equations  of  elasticity  (plane  strain)  in  the 
Cartesian  system  (x,y)  can  be  written  as: 

Lx(u,v)  «  0  (4.5) 


Ly(u,v)  ■  0 


(4.6) 


where 


Lx(u,v) 


(X+G>  & 


3u  3v 
3x  3y 


)  +  G  ( 


Ty 


3v 


)  +  G  ( 


,2  .2 

3  v  ,  3  v 

— ?  — 2 

3x^  dy* 


(4.7) 


Ly (u,v) 


(4.8) 
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with  X  and  G  being  Lame's  constants.  We  now  multiply 
equation  (4.5)  by  $  and  integrate  over  the  domain.  Similarly, 
we  multiply  equation  (4.6)  by  \p,  integrate  over  the  domain 
and  add  the  resulting  equations  to  obtain: 


5J[Vu,v 


)  ♦  +  L  (u,  v)  i|>]dA 


(4.9) 


Integrating  twice  by  parts,  as  explained  in  Chapter  2,  we 
obtain: 


[L  ((j),4))u  +  L  (<j>,^)v]dA-  \  [T  ($,^)u  +  T  (((),^)v]ds 

x  y  i  x  y 


[T  (u, v)  $  +  T  (u,v)i^]ds  =  0 
x  y 


(4.10) 


where : 


Vu'v)  *  <  It +  If  >  nx  * G(  H  "x +  W  V 


+  G(  n  +  |2.  n  ) 
3x  x  9x  y 


(4.11) 


Ty(u,v)  -  <&+§>  ny  +  G(  nx  +  ny) 


.  3u  _  .  3v  » 

+  Gl  t —  n  +  t —  n  ; 

3y  x  3y  y 


(4.12) 


are  the  tractions  corresponding  to  displacements  (u,v) 

and  n*(n  , n  )  is  the  outward  normal  to  the  surface.  In 
x  y 
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the  case  where  only  tractions  are  specified  on  the  boundary 
the  area  integral  in  (4.10)  vanishes,  since  the  extraction 
functions  satisfy  equilibrium  everywhere,  and  equation 
(4.10)  becomes: 


[T  (4>,i|0u  +  T  (<j>,^)v]ds  + 
A  y 


[Tx(u,v) 


+  Ty(u,v)^]ds  *  0  (4.13) 

The  integration  will  be  performed  separately  on  the 
circular  arc  and  the  rest  of  the  boundary  r.  On  the 
part  of  the  boundary,  where  tractions  are  specified, 

Tx(u,v)  and  Ty(u,v)  are  known: 

Tx(u,v)  -  X  (4.14) 


Ty(u,v)  «  Y  (4.15) 

On  the  part  of  the  boundary,  were  displacements  are  speci¬ 
fied,  <|)  and  t|>  vanish.  The  second  integral  in  (4.13) 
evaluated  around  the  contour  r  can  then  be  written: 


J 


[Tx(u,v)$  +  Ty  (u,v)i|>]ds  - 


J 


(X4>  +  Yip)  ds 


(4.16) 


This  integral  can  be  evaluated  numerically  from  the  data  of 
the  problem.  The  integration  around  the  circular  arc 
yields  an  expression  of  the  required  stress  intensity 
factor.  On  this  arc  T£  the  displacements  u  and  v  are  given 
by  equations  (4.1)  and  (4.2).  Their  derivatives  will  be  of 
the  form: 

=  Kx  r  1  (1^(6)  +  K2  r  2  02x(6)  +  ...  (4.17) 


fj)  *  *1.  r<1  uiy<e>  +  K2  D2yt9)  +  •• 


&  *  Ki  r<1'1  Vlx<9>  +  V2’1  V2x<0>  +  "• 


#  -  K1  r"X  1  Vly(9)  +  K2  ^  1  V2y(9)  +  '  \ 


(4.18) 

(4.19) 

(4.20) 


In  view  of  (4.11)  and  (4.12)  and  the  fact  that  on  the 
circular  arc  r£  the  outward  normal  n  has  components 

nx  *  -cos  9  (4.21) 


n^  ■  -sin  9  (4.22) 

the  tractions  T  (u, v)  and  T  (u,v)  can  be  written  in  the 

x  y 

form: 

K.-l  <  2~1 

T  (u,v)  »  Kl  r  1  Txl  (e)  +  K2  r  Tx2  (9)  +  . . . (4.23) 


”  »  *  •  *  »  *  *  *>"* 


•.  *»  % 


*.  *.  \ 
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K.-l  <0-^ 

T  (u, v)  *  K1  r  T  ,(0)  +  K2  r  2  T  2 (6)  +  ... (4.24) 


By  noting  that  on  r 


ds  *  r  d0 


(4.25) 


we  can  write  the  second  integral  in  (4.13)  around  F£  in 
the  form: 


J[Tx(u,v)*  +  Ty(u,v)t^]ds  ®  +  k2  M2  (r) 


+  k3  M3(r)  +  K4  M4(r) 


(4.26) 


where 


Mi  -  ( “ 
J  -a 


[Txl(0)  «(6)  +  T  1(0)0(6)]de 


(4.27) 


M2  (r) 


-<  ^+tc  2  f  a 

r  \ 

J  -a 


[tx2(9)  0(0)  +  Ty2(0)  Y(0)]d0 


(4.28) 


M3(r)  «  r 


■S: 


[Txl(0)0b  +  Tyl(0)0b]d0 


(4.29) 


M4(r)  -  r 


tTx2(e)0b  +  Ty2{e)^de  {4'30) 


The  higher  order  terms  are  disregarded  in  what  follows , 
since  two  terms  are  sufficient  for  the  discussion  of  all 
the  relevant  issues.  We  now  observe  that  is  independent 
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of  the  radius  r.  This  is  the  reason  that  the  extraction 
functions  are  chosen  as  the  eigenfunctions  corresponding 
to  the  negative  eigenvalue  -k^.  Otherwise  would  also 
depend  on  r.  It  can  be  seen  that  M3 (r)  and  M4(r)  vanish 
as  r  approaches  zero . 

In  view  of  (4.3)  and  (4.4)  the  derivatives  of  the 
extraction  functions  can  be  written  in  the  form: 


9$  _ 
3x 


$  (6)  +  (J) 

x 


bx 


*  (6)  +  <J> 
Yx(0)  +  ^ 


by 


bx 


(4.31) 

(4.32) 

(4.33) 


|f  -  r  Kl  1  »y(9)  +  *by  U-34) 


where  $b2,  $by,  ^bx,  ^by  are  again  smooth  functions.  The 
tractions  corresponding  to  displacements  4>  and  \p  can  be 
written  in  the  form: 

-<  -i 

TX(<P,V)  -  r  X^(9)  +  Xb  (4.3 

-<,-1 

Ty  (<P  ,ip)  -  r  Y^(6)  +  Yfa  (4.3 

with  X.  and  Y.  being  again  smooth  functions.  The  first 
D  o 

integral  around  T£  in  (4.13)  can  then  be  written  as: 


rt.ViVA.V, 


.VaWuV: 
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We  observe  that  expression  Mg  is  also  indpenedent  of  r. 

The  choice  of  the  extractions  functions  <t>  and  ip  as  the 
eigenfunctions  corresponding  to  the  negative  eigenvalue 
makes  both  Mg  and  Mg  independent  of  r.  Since  and 
are  smooth  functions,  expressions  M7 (r)  and  Mg (r)  vanish  in 
the  limit  as  r  approaches  zero.  Expressions  M2 (r)  and 
Mg(rl  require  further  consideration.  In  the  case  where 
is  the  first  eigenvalue,  that  is  <  <2»  then  the  quantity 
(-<^+<2)  is  positive  and  both  M2(r)  and  Mg (r)  vanish  as  r 


tends  to  zero.  However,  if  a  higher  order  stress  intensity 
factor  is  extracted,  say  K2,  the  exponent  of  r  in  expres¬ 
sions  M2(r)  and  Mg (r)  will  be  negative.  These  expressions 
then  become  unbounded  as  r  tends  to  zero,  unless  the 
integrals  with  respect  to  0  vanish.  It  will  now  be  shown 
that  they  indeed  vanish. 

Consider  the  body  shown  in  figure  (4.4),  which  is 
bounded  by  the  straight  segments  r 2>  r4  on  the  faces  of 
the  notch  angle  and  the  circular  arcs  T ^ of  epicentral 
angle  2a  centered  at  the  notch  tip.  Applying  equation 
(4.13)  for  this  body  with  u,  v  as  the  second  terms  of  (4.1) 
and  (4.2)  respectively  and  <p ,  ip  as  the  singular  terms  of 
(4.3)  and  (4.4)  we  obtain: 


r  ^<1+k2-i 


Cx„<9)  n2(e)  +  v2(e)  ]rae 


ri+r3 


f  —k.+Kj—I 

\  r  [Tx2(e)(He)  +  Ty2(e)'Me)  ]rde  =  0 

(4.42) 


r  +r 
1  3 


Both  integrals  vanish  on  the  straight  segments  ^ ,  r4 
because  the  traction  components  Tx2»  T  2  corresPon<*ing  to 
the  eigenfunctions,  as  well  as  the  tractions  components 
X  ,  corresponding  to  the  extraction  functions  vanish 

on  the  faces  of  the  angle. 
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Subscripts  I  and  II  refer  to  modes  I  and  II  respectively. 
The  stresses  in  the  Cartesian  system  (x,  y)  can  be 
expressed  through  the  transformation  formulas: 


2  2 

a  st  a  cos  6  +  ac  sin  0-2  t  D  sinQ  cos0  (5.13) 

x  r  w  ro 


2  2 

a  =  a  sin  0  +  aD  cos  0+2  t  a  sin©  cos0  (5.14) 

y  r  o  ry 


2  2 

T  *  (a  -aa)sin0  cos  ©  +  t  Q (cos  0-sin  0)  (5.15) 

xy  r  y  ru 


The  edges  of  the  notch  are  free  from  tractions  with 

stresses  a.  and  t  q  being  zero  there.  On  any  other  plane 
u  X  u 

through  the  body  with  an  outward  normal  n  *  (n  ,  n  )  the 

x  y 

tractions  are  given  in  terms  of  stresses  through  the 
transformation  law: 


X 


Y 


ax  cosa  + 


sina  + 


sina 


cosa 


(5.16) 

(5.17) 


An  alternative  way  of  obtaining  expressions  for  the 
tractions  is  to  obtain  first  the  derivatives  of  the 
displacements  and  then  express  the  tractions  in  terms  of 
these  derivatives  as  follows: 
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1/2  * t  * 

CjCZk)  Trei  =  Kj  r  <lC  (<I+D  sin(KI+l)6  + 


<T-1 


+  Cj  (<I»1)  sin^-l)  6] 


(5.7) 


1/2  ^  tt ~ 

Ctt  ( 2tt  )  a_TT  =  -K,,  r  [K  n  n+D  sin  (<  j  j+D  9  + 


II 


rll  *'II 


+  CII<H  (<H-3)  sin(<n-D  8]  (5.8) 


1/2  tt"^ 

Cn(27r)  aen  *  Kn  r  <II  (|CiI+l)  C  sin  (k  ii+l)  9  + 


eii  “  ii 


+  CIIsin(<II-l)  e] 


(5.9) 


1/2  ^  tt“1 

cn(27r)  Tren  *  “Ku  r  (|CII+i)c°s(|CII+i)e  + 


+  Cjj  (Kjj-1)  cos  (iCjj-I)  e] 


(5.10) 


where 


cos (<j+l) a 


cos  (Kj-1)  a 


(<2+1)  sin(<I+l)  a 


(<^-1)  sinU^-l)  a 


(5.11) 


sin  (< 2j+l)  o 


'II 


sinU^-l)  a 


(<22+1)  cos  (k  +1)  a 


(Kjj-l)  cos  [k~^TT  a 


(5.12) 


■  v 

j 


C*  •  *  •  *  •  *  *  . 
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2GCi(2t7)1/2  ux  -  Kx  r  I[  (C;[B-<I-1)  cosk^  - 

-  2kiCicos0cos (Kj-l) 6]  (5.1) 

2GCj  ( 2  tt  ) 1/2  vx  -  Kx  r  I[  (CjB+Kj+1)  sinic^  - 

-  2<ICIsin0cos (Kj-l) 0]  (5.2) 


2GCI;E(27r) 


1/2 


u 


II 


2GCjI  ( 2ir ) 


1/2 


II 


CI(2ir) 


1/2 


CjCir)1/2  n9l  - 


-  Kjj  r  I3:[  (CjjD-Kjj-DsinKjje  + 

+  2icIICIIsin0cos  (Kjj-l)  0]  (5.3) 

^  tt 

*  KXI  r  j"l[-(CiiD+<ii+1)cosicii0  - 

-  2<IjCIIcos0cos (Kjj-1) 0]  (5.4) 

<  -1 

-Kj  r  [<I  (<I+1)  cos  (<x+l)  0  + 

+  CjKj  (Kj-3)  cos  (Kj-1)  9]  (5.5) 

<  T  “ 1 

KI  r  1  <!  (<!+l)  Ccos  (<I+D  6  + 

+  CjCOS (<x-l) 9]  (5.6) 


5.  THE  MODEL  PROBLEM 


The  numerical  performance  of  the  extraction  method 
will  now  be  established  through  numerical  experiments 
based  on  a  model  problem.  The  model  problem  contains  a 
reentrant  corner  and  a  solution  is  selected  which  satis¬ 
fies  (1)  the  Navier  equations  and  (2)  the  boundary 
conditions  at  the  reentrant  corner.  On  the  other 
boundaries  the  tractions  corresponding  to  the  exact  solution 
are  specified.  In  this  way  the  exact  solution  has  the 
main  characteristics  of  typical  practical  problems  (i.e. 
ones  with  reentrant  corners) .  It  is  therefore  suitable 
for  benchmark  studies.  No  problem  with  reentrant  corners 

and  possessing  an  exact  solution  has  been  reported  in  the 

» 

literature.  In  all  the  reported  solutions  a  numerical 
approximation  is  made  at  some  stage,  in  which  the  error 
is  unknown.  Thus  precise  convergence  studies  cannot  be 
performed . 

Let  us  consider  the  two-dimensional  domain  of 
infinite  dimensions  with  a  notch  of  total  solid  angle 
equal  to  2a.  As.  it  was  explained  in  Chapter  3  the  stress 
and  deformation  fields  are  given  by  series  expansions  in 
which  the  typical  terms  are  as  follows: 


The  same  will  happen  when  more  terms  are  kept  in  the 
expansion.  Equation  (4.76)  then  reduces  to: 


By  choosing  in  turn  the  extraction  functions  to  be  the 
eigenfunctions  with  the  negative  eigenvalue  -<1  and 
corresponding  to  the  eigenfunctions  that  have  as  stress 
intensity  factors  the  constants  K2,  ,  K4  respectively 

in  equations  (4.54),  (4.55),  we  obtain  three  more  equations 
of  the  same  form.  We  have  then  a  system  of  four  equations 
in  four  unknowns  to  determine  K^,  K2,  ,  and  K4 . 
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where : 


Mi 


4 

I 

i«l 


hij (0)de 


(4.77) 


(4.78) 


We  have  shown  that  the  integral  in  (4.78)  must  be  indepen¬ 
dent  of  r.  Since  the  functions  c(r),  s(r),  c^r),  s^r) 
are  linearly  independent,  this  can  happen  only  if: 

M2  *  M3  »  M4  «  M5  «  0  (4.79) 

So  far  only  one  term  in  the  expansion  has  been  considered. 
The  inclusion  of  further  terms  will  give  rise  to  an 
expression  of  the  form: 

+  c(r)M2  +  s(r)M3  +  c1(r)M4  +  s1(r)M5 

K  <  1C  < 

+  r  3c(r)Mg  +  r  3s(r)M7  +  r  3c1(r)Mg  +  r  3s1(r)M9 

(4.80) 

as  the  right-hand- side  of  (4.76).  It  can  be  shown  that 

K  3  K  3  IC  3 

the  new  functions  of  r  (c,  sf  r  s,  r  ci'  r  s^) 

sure  also  linearly  independent.  Thus  expression  (4.79)  can 
be  independent  of  r  only  if : 


M.  - 


1 


0  i  «  2 


9  •••  9 


9 


(4.81) 


Fi(r,9)  =  9li(0)  +  c(r)  g2i(8)  +  s(r)  g3i(0) 


+  c1(r)  g4i ( 6 )  +  s1(r)  g5i(8)  (4.73 

where  now  c(r),  s(r),  c^(r),  s^(r)  are  functions  of  r  that 

are  linearly  independent.  The  other  three  products 

T  (u,v)i|/ds,  T  (<j>,i|/)xds  and  T  (<|>,^)vds  are  of  the  same 

y  x  y 

form  and  we  can  write  for  the  contour  r^: 

“  [Tx(<J>,i|>)u  +  Ty (4>,i|>,v]ds  +  [Tx(u,v)<()  +  Ty (u,v)^]ds  * 

4 

-  Z  K.  H. (r, 6) d8  (4.74 

i»l  1  1 

where : 

Hi(r,8)  »  hli(8)  +  c(r)h2i(8)  +  s(r)h3i(8)  + 

+  c1(r)h4i(8)  +  s1(r)h5i(8)  (4.75! 


Integrating : 


-i 


[Tx(<J>,i|/)u+Ty(<t>,'}>)  v]ds  + 


)CVn' 


v)  4>+Ty  (u,v)  ^]ds 


e  e 

M,  +  c(r)M-  +  s(r)M,  +  c, (r)M.  +  s,  (r)M-  (4.761 


T  (u,v)4>ds  “IK.  F.  (r,0)d0 
x  i«l  1  1 


where 


(4.66) 


Fi(r,0)  »  Gn(9)  +  c(r)G2i(0)  +  s(r)G3i(0)  + 

+  c2(r)G4i(0)  +  s2(r)G5i(0)  +  c  (r)  s  (r)  Ggi  (0) 

(4.67) 


By  denoting; 

c^(r)  “  cos(2<2^nr)  (4.68) 

s1(r)  «  sin(2ic2Jtnr)  (4.69) 

we  can  write; 

c2(r)  «  cos2(<2*'nr)  *  \  [c^(r)-l]  (4.70) 

s2(r)  »  sin2(<2^i^)  *  j  [c^(r)+l]  (4.71) 

c(r)s(r)  ■  cos  (K^lnr)  sin(K2inr)  “  J  s^r)  (4.72) 


and  rearrange  (4.67)  to  obtain; 


w 
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<,-l 


Tx(u,v)  -  r  1  K1CTxll(e)+c(r)Txl2 (9)+s(r)Txl3 (0) ]  + 


■c,-l 


+  r  1  K2[Tx21(6)+c(r)Tx22(e)+s(r)Tx23(9)] 


<,-l 


+  r  1  K3[Tx3^ (0 ) +c (r) Tx32 (9 ) +s (r) Tx33 (0 )  ]  + 


K,-l 


+  r  1  K4[Tx41(0)+c(r)Tx42(0)+s(r)Tx43(0)]  (4.62) 


k,-1 


(u,x)  »  r  1  K1[T  4 (0)+c (r)T  12 (0)+s (r)T  13  (0) ]  + 


<l“1 


+  r  x  K2[T  (6)+c(r)T  22(0)+s(r)T  23(0)] 


<•.-1 


+  r  1  K3[T  3l(0)+c(r)T  32(0)+s(r)T  33(0)]  + 


<,-l 

+  r  x  K4[Ty41(0)+c(r)Ty42(0)+s(r)Ty43(0)]  (4.63) 


-<.-l 


Tx  (4>,^)  ■  r  1  Exx  (8)  +  c(r)  X2(0)  +  s(r)  X3  (0)3  (4.64) 


-<,-l 


T  (<fr,i|»)  -  r  1  [Y1(9)  +  c  (r)  Y2(6)  +  s(r)  Y3(0)](4.65) 


We  are  now  in  a  position  to  evaluate  the  integrals 
appearing  in  (4.13)  around  the  circular  arc  r  .  The 
product  Tx(u,v)<))ds  can  be  written  as: 


+  r 


1 


K3[v3i(e)  +  c(r)  V32  (9)  +  s(r)  V33(0)]  + 

+  r  1  K4[V41(6)  +  c (r)  V42(0)  +  s (r)  V43(0)] 

(4.55) 


where : 


c(r)  =  cos  («2&nr)  (4.56) 

s(r)  «  sin  (<22.nr)  (4.57) 

The  corresponding  extraction  functions  for  extracting  <1 
are: 

<j>  »  r  i[4>1(9)  +  c(r)  <t>2(0)  +  (4.58) 

n»  -  r  1[H'1(0)  +  c  (r)  ¥2(0)  +  s  (r )  V3  (0)  ]  (4.59) 


Since 

— [cos(<2inr)]  ■  -  <2  r"1  sin(<2Jlnr)  (4.60) 

[sin (<2lnr)  ]  *  <2  r”1  cos(<2inr)  (4.61) 


the  corresponding  tractions  can  be  written  in  the  form: 
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As  it  was  explained  in  the  previous  chapter,  to  the  posi¬ 
tive  real  part  there  correspond  four  independent  eigen¬ 
functions,  each  with  its  own  stress  intensity  factor  and  a 

S 

singularity  of  the  order  r  .  The  smooth  part  also 
contains  terms  of  the  form: 

sin  (<2Jlnr)  ,  cos(<25,nr)  (4.53) 

which  means  that  the  eigenfunctions  can  be  written  in 
the  form: 

£  l 

u  -  r  x  K1[U11 (6)  +  c  <r)  O12(0)  +  s(r)  U  3  (0) ]  + 

k  -l 

+  r  K2[021(6)  +  c (r)  U22(0)  +  s(r)  023 (0) ] 

+  r  K3[U31(0)  +  c (r)  U32(0)  +  s(r)  U33(0)  + 

+  r  K4[U41(0)  +  c  (r)  U42(0)  +  s(r)  O43(0)] 

(4.54) 

v  =  r  K1[V11(0)  +  c  (r)  V12(9)  +  s  (r)  V13(0)]  + 

+  r  K2[V21(0)  +  c  (r)  V22(9)  +  s(r)  V23(9)] 


Equation  (4.13)  then  becomes: 


Ki(-M2+M5)  *  \  x<J>ds  +  1  Y  ds  -  1  [Tx(4>,^)u  + 


+  Ty  (4>,<{>)v]ds 


(4.50 


When  displacements  cure  prescribed  on  some  part  of  the 
boundary,  the  area  integral  is  also  present  and  (4.50) 


becomes : 


^(-Mj+Mg)  »  ^  a^[Lx(4>,^)x  +  LyU,i|>)v]dA  + 


+  \  X<j)ds  +  \  YiJ;ds  -  \  [Tx(<j>,(j;)x  +  Ty  (<J>  ,ip)  v]ds 


(4.51 


In  order  to  obtain  an  approximation  to  we  substitute  in 
(4.50)  or  (4.51)  u  «  UpE  and  v  «  vpE. 

The  case  of  complex  eigenvalues : 


,der  the  case  of  a  pair  of  complex  conjugate 
e  form: 


+  i  k. 


(4.52) 


L. 
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P 


I 


or: 

-M2(r3)  +  Mg(r3)  -  -M2(r3*)  +  Mg(r3*)  (4.45) 

* 

where  r3  is  the  radius  of  the  arc  r3  and  r3  is  the  radius 
* 

of  the  arc  r3  .  This  means  that  the  expression  (-M2+Mg ) 
takes  the  same  value  for  any  value  of  the  radius  r.  By 
writing  this  expression  in  the  form: 


-<.+<_  fa 

-M2  +  M6  *  r  \ 

J  -a 


“Kl+K2 

F ( 0 ) d0  *  r  x  z  C 


(4.46) 


it  is  obvious  that  this  can  be  independent  of  r  only  if 
the  integral  is  zero  (the  eigenvalues  <2  are  distinct) , 
therefore : 


-M2  +  Mg  *  0  (4.47) 

When  more  terms  are  kept  in  the  asymptotic  expansion  the 
corresponding  expression  will  be: 

-K.+K-  -<.+<_  — < .  +IC  . 

-M2  +Mg»r  x  xC2+r  1  JC3+r  1  C4  +  . . . 

(4.48) 

where  the  eigenvalues  «2,  «3,  <4  are  all  different  from  <^. 
This  expression  can  be  independent  of  r  only  if: 

C2  *  C3  *  C4  ■  0  (4.49) 

consequently,  we  obtain  again  equation  (4.47). 


L  •  •  *  • 


AV.V. 


2cnG<2*>  '  if  -  KII  r  kiiL(cii-,cii-1) 

sin(<II~l)9  +  2  CjjCKjj-I) 
sinScos (<II"2) 8]  (5.22) 

•i  ,,  3U  <TT-1 

2CIIG(2it)  3y  “  KII  r  <II^  {CIID“<Il“1) 

cosUjj-l)©  -  2<IICIIsin9 

sin (<II"2)  8  +  2CI1cos8cos(icII-2)  0 

(5.23) 

.  3v__  <tt_1 

2ClIG( 2ir>l/2  -  KXI  r  '„t-(Cn”«„+l) 

•  COs(kii-1)8  -  2< ^ jCj jCos8 
cos (<II"2) 8  +  2CIIsin8sin (<II-2) 8 


(5.24) 


-123- 


i  /o  ^tT  ^ 

2CIIG(2,'>  -  KjI  r  «n[(Cn^n*l) 


sinte^-DS  +  2cu  (|CII"D  cosQ 


sin (k^-2)  0] 


(5.25) 


XI  - 


(X+2G) 


Yj  «  (X+2G) 


3u 

3v 

3u 

3v 

T3 r  nx  + 

X 

3y 

nx  +  6( 

~  + 

I 

t r 

3v 

3u 

3u_ 

3v 

Ty“  ny  + 

X 

X 

3x 

ny  +  G( 

TF  + 

I 

3x 

)n 


)n. 


(5.26) 


(5.27) 


kII 


3u__  3v__  3u__  3v__ 

‘i+2G>  -#  nx  +  x  -#  “x  +  G<  -#  +  -S?  )ny 

(5.28) 


3v 


YII  *  (X+2G> 


II 


\  3Ul1  n  i  rf  3UjI  l  3V]CI  In 

w  y  “Sx"  y  G(  “5x"  )nx 

(5.29) 


Let  us  now  cut  out  from  the  infinite  domain  a  body 
of  finite  dimensions  which  contains  the  notch ,  but  is 
otherwise  arbitrary.  If  we  apply  to  this  body  the 
tractions  specified  by  (5.26)  to  (5.29),  then  the  state 
of  stress  and  deformation  will  be  given  by  (5.1)  to 
(5.10)  in  accordance  with  the  uniqueness  theorem  of 
elasticity.  Since  the  body  was  in  equilibrium  when 
embedded  in  the  infinite  domain,  it  must  be  in  equilibrium 
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under  the  tractions  on  the  cut-out  surface.  When  these 
tractions  are  then  applied  as  external  loading,  this 
loading  is  self-equilibrating.  The  displacement  boundary 
conditions  that  have  to  be  enforced  are  then  the  ones 
necessary  to  prevent  rigid  body  motions .  In  the  problems 
that  we  analyzed  we  chose  to  fix  both  displacements ' at 
the  notch  tip  and  the  u  displacement  at  the  upper  right 
corner . 

We  can  then  predetermine  the  singular  elastic  field 
that  exists  throughout  the  body  by  arbitrarily  choosing 
the  generalized  stress-intensity  factors.  The  applied 
tractions  necessary  to  induce  this  field  are  then  given 
by  (5,26)  to  (5.29).  Any  number  of  terms  can  be  retained 
in  the  expansion.  Thus  we  can  create  problems  for  which 
the  exact  values  of  the  generalized  stress  intensity 
factors  of  any  order  are  known.  In  the  problems  studied 
herein  only  the  first  terms  in  the  expansions  were  kept 
and  the  stress  intensity  factors  and  were  chosen 
to  have  an  exact  value  of  unity.  The  combined  loads 
corresponding  to  mode  I  and  mode  II  were  applied  in  each 
case.  Three  problems  were  analyzed  with  a  total  solid 
angle  of  360,  330  and  270  degrees  respectively  as  shown 
in  figures  5,1,  5.2,  5.3.  The  last  case  is  one  that 
occurs  very  frequently  in  practice  representing  a  right 
angle  cut-out  (L-shaped  domain) .  The  results  are  listed 


-128- 


in  tables  5.1,  5.2,  5.3  and  are  shown  graphically  in 
figures  5.4  through  5.12.  The  theoretically  predicted 
rates  of  convergence  are  also  shown  in  the  figures.  The 
rate  of  convergence  in  energy  is  equal  to  twice  the 
smallest  eigenvalue,  which  is  the  first  mode  I  eigenvalue. 
For  the  stress  intensity  factors  the  theoretical  rate  of 
convergence  is  equal  to  twice  the  corresponding  eigenvalue, 
since  the  tractions  corresponding  to  mode  II  have  no 
influence  on  the  mode  I  stress  intensity  factor  and  vice 
versa.  It  can  be  seen  from  the  figures  that  the  theore¬ 
tical  rates  of  convergence  are  verified  numerically  and 
even  for  the  coarse  mesh  chosen  a  relatively  small  number 
of  degrees  of  freedom  is  required  to  give  answers  within 
required  engineering  accuracy  (usually  1  to  5  percent) . 


Convergence  results  for  the  model  problem  of  figure  5.2  (solid  angle  =  330°) 


Convergence  results  for  the  model  problem  of  figure  5.3  (solid  angle  =  270°) 


NUMBER  OF  DEGREES  OF  FREEDOM 


Figure  5.4 

Energy  convergence  for  model  problem  of  Figure  5 . 1 

(solid  angle  =*  360°) 


REL.  ERROR  IN  ENERGY  (PERCENT) 


RELATIVE  ERROR  (PERCENT) 


■('i  ■■  ■  «  my  w\ i  ■;  ^  i  ^  '^v. I1  ■  v  ui j  i »  ■■  vp.'mt  * 
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Figure  5.7 

Convergence  of  Kj  for  model  problem  of  Figure  5.1 
(solid  angle  =  360°) 


1000 
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vicinity  of  reentrant  corners  and  therefore  can  be 
utilized  in  research  in  the  area  of  crack  initiation. 

Direct  methods  for  stress  computation,  whether  based  on 
the  finite  element  method  of  other  numerical  methods, 
have  the  disadvantage  that  in  the  vicinity  of  singular 
points  they  incur  large  errors,  therefore  they  are  not 
well  suited  for  correlation  between  experimental  and 
analytical  data.  In  order  to  realize  the  full  benefits 
of  the  extraction  method,  extensive  experimentation  must 
be  performed  with  the  goal  to  determine  relationships 
between  elastic  stress  field  parameters  and  failure 
initiation  events.  Preliminary  experiments  with  brittle 
epoxy  resin  specimens  were  conducted.  Monotonicity  was 
found  to  exist  between  the  stress  intensity  factor  and 
failure  initiation  for  a  wide  range  of  solid  angles.  Once 
theories  of  failure  initiation,  based  on  elastic  stress 
field  parameters,  are  developed  and  proven  the  extration 
method  presented  herein  will  provide  a  basis  for  the  design, 
analysis  and  certification  of  structures  with  geometric 
singularities. 


» 
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7.  SUMMARY  AND  CONCLUSIONS 

The  generalized  influence  function  method  for  the 
extraction  of  stress  intensity  factors  in  plane  elasti¬ 
city  was  evaluated  from  the  point  of  view  of  its  potential 
utility  in  engineering  applications.  The  method  was 
implemented  and  applied  in  conjunction  with  the  p-version 
of  the  finite  element  method.  The  method  yields  the 
amplitudes  of  all  terms  in  the  asymptotic  expansion  of 
the  linear  elastic  field  in  the  neighborhood  of  sharp 
corners  and/or  sudden  changes  in  boundary  conditions . 

In  this  investigation  only  reentrant  corners  with  free- 
free  boundary  conditions  were  considered. 

A  model  problem  for  which  the  exact  solution  is 
known  was  used  to  test  the  implementation  and  for 
evaluating  the  method.  Convergence  studies  were  performed 
and  the  theoretically  predicted  rates  of  convergence  were 
verified  numerically.  The  method  was  found  to  be  accurate 
and  reliable.  Even  with  coarse  meshes  and  relatively  small 
number  of  degrees  of  freedom,  it  was  possible  to  compute 
the  stress  intensity  factors  to  within  a  few  percent  rela¬ 
tive  error. 

The  engineering  importance  of  the  extraction  method 
presented  herein  is  that  it  provides  good  qualitative 
and  quantitative  understanding  of  the  stress  field  in  the 
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Table  6 . 1 

Results  of  experiments 


Solid 

Anglo 

First 

Eigenvalue 

(X) 

Specimen 

Specimen  dimensions  (bd) 

Failure 
Load  (N) 

Computed 

(MPa  m1*X) 

Width 

Ligament 

Thickness 

360* 

m 

2.5 

38.1 

28.4 

3.2 

657 

0.805 

m 

38.1 

28.6 

3.3 

804 

0.849 

38.1 

28.6 

3.3 

824 

0.882 

330* 

0.5015 

19 

38.1 

28.6 

3.3 

785 

0.839 

4+ 

31.3 

22.2 

3.2 

647 

0.863 

- 

m 

38.1 

28.7 

3.3 

1491 

1.645 

0.5122 

wm 

38.1 

28.7 

3.3 

1549 

1.724 

m 

38.1 

28.7 

3:3 

1461 

1.614 

i 

■■ 

3.1 

38.1 

26.7 

3.3 

1451 

1.955 

270* 

3.2 

38.1 

26.7 

3.3 

1383 

1.841 

Mm 

3.3 

38.1 

26.7 

3.3 

1461 

1.968 

IS 

38.1 

29.0 

3.3 

1648 

2.602 

240* 

0.6157 

B9 

38.1 

29.0 

3.3 

1608 

2.530 

38.1 

29.0 

3.3 

1589 

2.509 

wm 

38.1 

28.2 

3.3 

3452 

9.073 

KB 

38.1 

28.2 

3.3 

3364 

8.800 

210* 

0.7520 

mm 

38.1 

28.2 

3.3 

3609 

9.463 

m 

38.1 

28.2 

3.2 

3383 

9.315 

Figure  6.2 

Typical  test  specimen 


Eighteen  speciments  were  tested,  a  typical  specimen 

is  shown  in  figure  6.2.  The  dimensions  of  the  specimens, 

the  failure  load  and  the  stress  intensity  factors  at 
* 

failure  (K^)  are  given  in  Table  6.1.  The  dimension  of 
* 

Kj  depends  on  the  strength  of  the  singularity,  which  is 

the  corresponding  eigenvalue.  Due  to  the  brittleness  of 

the  material,  the  failure  initiation  and  failure  events 

occured  at  virtually  the  same  load  level.  The  experimental 

results  are  plotted  in  figure  6.3.  Monotonicity  of  the 

critical  value  of  the  generalized  stress  intensity  factor 
* 

Kr  is  observed.  It  is  hoped  that  future  experimental  work 
will  establish  precise  failure  initiation  criteria  in 
terms  of  generalized  stress  intensity  factors  and  clarify 
their  range  of  applicability. 
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the  range  of  solid  angles  within  which  the  hypothesis  is 
valid.  In  the  context  "valid"  is  meant  in  the  sense: 
"sufficiently  accurate  for  the  purpose  of  predicting 
crack  initiation  in  practical  applications”. 

This  is  one  of  several  possible  hypotheses.  An 
extensive  amount  of  experimental  work  is  required  in  order 
to  identify  the  minimum  number  (and  kind)  of  material 
parameters  that  need  to  be  measured  in  order  to  be  able 
to  predict  crack  initiation  events  with  about  the  same 
accuracy  and  reliability  that  crack  extension  is  pre¬ 
dictable  today.  The  ability  to  compute  with  high  accuracy 
the  elastic  stress  field  parameters  for  any  sharp  notch 
configuration  allows  alternate  hypotheses  to  be  made  and 
makes  this  a  very  promising  research  area. 

The  material  used  for  the  experiments  was  a  two 
component  epoxy  system:  83  percent  Shell  Epon  828  resin, 
17  percent  Shell  Epon  curing  agent  Z.  Modulus  of  elasti¬ 
city  of  this  material  was  4.3  GPa  and  Poisson's  ratio 
0.3P-0.39  (based  on  ultrasonic  measurements).  The  resin 
was  cured  24  hours  at  room  temperature  and  12  hours  at 
65  degrees  Celsius.  The  equipment  was  an  INSTRON  testing 
machine  with  a  500  kgf  load  cell.  The  test  arrangement 
is  shown  in  figure  6.1. 
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6 .  PRELIMINARY  EXPERIMENTS 


A  number  of  preliminary  experiments  were  performed 
in  order  to  obtain  baseline  data  for  the  hypothesis  that, 
in  the  case  of  sharp  reentrant  corners,  crack  initiation 
is  controlled  by  the  intensity  of  the  stress  singular 
term. 

This  hypothesis  implies  that  crack  extension  and 
crack  initiation  are  not  fundamentally  different  phenomena; 
both  events  can  be  correlated  with  the  intensity  of  the 
stress  singular  term  at  the  reentrant  corner,  provided 
that  the  solid  angle  is  not  much  smaller  than  360  degrees. 
Of  course,  it  is  very  likely  that,  in  the  case  of  smaller 
solid  angles,  the  intensities  of  the  higher  order  terms 
may  also  have  to  be  considered. 

Assuming  this  hypothesis  to  be  valid,  one  would 
expect  the  intensity  of  the  stress  singular  term  at 
failure  initiation  to  be  some  monotonic  function  of  the 
solid  angle.  Of  course,  monotonicity  itself  is  not 
sufficient  evidence  that  the  hypothesis  is  valid.  A  number 
of  carefully  designed  experiments  involving  various  stress 
fields  must  be  performed  and  the  outcome  (i.e.  the  crack 
initiation  event)  successfully  predicted  from  baseline 
data,  representing  critical  material  stress  intensities, 
in  order  to  establish  the  validity  of  the  hypothesis  and 
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